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ABSTRACT 

We recover the joint and individual space density and surface brightness distribution(s) of 
galaxies from the Millennium Galaxy Catalogue. The MGC is a local survey spanning 
30.9 deg 2 and probing approximately one-two mag arcsec' 2 deeper than either the Two- 
Degree Field Galaxy Redshift Survey (2dFGRS) or the Sloan Digital Sky Survey (SDSS). 
The MGC contains 10 095 galaxies to B Mac < 20 mag with 96 per cent spectroscopic com¬ 
pleteness. For each galaxy we derive individual /f-corrections and seeing-corrected sizes. 
We implement a joint luminosity-surface brightness step-wise maximum likelihood method 
to recover the bivariate brightness distribution (BBD) inclusive of most selection effects. In¬ 
tegrating the BBD over surface brightness we recover the following Schechter function pa¬ 
rameters: 4>* = (0.0177 ± 0.0015)/i 3 Mpc -3 , M* - 5 log ft = (- 19.60 ± 0.04) mag 
and a = — 1.13 ± 0.02. Compared to the 2dFGRS dNorberg et al J2002t) we find a consistent 
M* value but a slightly flatter faint-end slope and a higher normalisation, resulting in a final 
luminosity density jbj = (1.99 ± 0.17) x 10 % hL^ Mp c~ 3 — marginally high er than, but 
consiste nt with, the earlier 2dF GRS dNorberti et all2002l) . ESP dZucca et all 19971) . and SDSS 
z = 0.1 (Blan ton et alJl2003ah results. The MG C is inconsistent with the SDSS z = 0 result 
(+3 <t) if one adopts the derived SDSS evolution dBlanton et all2003alh 

The MGC surface brightness distribution is a well bounded Gaussian at the M* point 
with <j>* = (3-5 ± 0.1) x 10 ~ 2 h 3 Mpc -3 , fi e * = (21.90 ± 0.01) mag arcsec -2 and 
din r c = 0.35 ± 0.01. The characteristic surface brightness for luminous systems is invari¬ 
ant to Mb ucc — 51og/i — —19 mag faintwards of which it moves to lower surface bright¬ 
ness. The surface brightness distribution also broadens (<7i n -0.5-0 .7) towards lower 
luminosities. The luminosity dependence of o\ u r, provides a new constraint for both the the¬ 
oretical development (Dalcanton, Spergel & Summers 1997; Mo, Mao & White 1998) and 
numerical simulations (e.gJColejfcLace’^1996|) which typically predi ct a mass-independent 
crin R e — 0-56 ± 0.04 (see I Vitvi'tskaetail20() 2 and Bullock et al 1200 111 . 

Higher resolution (FWHM -C 1 arcsec) and deeper (//, IM 26 mag arcsec -2 in the B- 
band) observations of the local universe are now essential to probe to lower luminosity and 
lower surface brightness levels. 

Key words: galaxies: general - galaxies: fundamental parameters - galaxies: luminosity 
function, mass function - galaxies: formation - galaxies: evolution - astronomical data bases: 
catalogues 


1 INTRODUCTION 

The recognised convention for representing the overall galaxy pop¬ 
ulation is the luminosity distribution or space density of galaxies 
(see review by Binggeli, Sandage & Tammann 1988 and more re- 


* E-mail: spd@mso.anu.edu.au 
f PPARC Visiting Fellow 


centIv iDriver & De Pronri J20()3fl . This distribution is typically de¬ 
rived from a magnitude limited redshift survey which is corrected 
for Malmquist bias (e.g., Efstathiou, Ellis & Peterson 1988) and fit¬ 
ted with the three parameter Schechter function ISchechteil 197(1 ). 
Generally, the Schechter function is found to be a formally good 
representation, alth ough some surveys (e.g., IZucca et*a ff 1997: 
\1 arz k e et al J1 1 9981 ) have hinted at an upturn at the faintest lim- 
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its, * 1 where the statistics_become poor (typically Mb — 5 log h w 
— 16 mag, see also lBlanton et alJI2005all . The process of fitting a 
luminosity function (LF) reduces the galaxy population to three 
crucial numbers: the characteristic luminosity, M* (or L* ), the nor¬ 
malisation, <j>*, and the faint-end slop e, a. 

Since the pioneering work of lDavis & Huchra ( 1982 1 numer¬ 
ous measurements of these three crucial parameters have been 
made (see lLiske et alJ2003i) . Among the most recent are those de¬ 
rived from the Two- Degree Field Galaxy Redshift Survey (2dF- 
GRS; (Norberg e t alj [2002}) and the Sloan Digital Sky Survey 
(SDSSl&2~ Blantone^ni200lt 2003a). However, although indi¬ 
vidual surveys yield Schechter function parameters to high ac¬ 
curacy, a comparison between independently published values 
shows excessiv ely large variations, indicative of strong sy stem- 
atic errors (see [Cross et akll200ll: Elske et alJl20oA Ide Lapparenl 
12003 1. Cross & Drived 12002^ discuss the impact of surface bright¬ 
ness selection effects on measurements of the LF (see also 
ISpravberrv et all 1 9971 : lDalcantonlll998t) and find that such effects 
may p rovide an explanation for the variation seen. In lLiske et alJ 
120031) we compared a number of recent B-band Schechter func¬ 
tion values and used the precision number counts of the Millen¬ 
nium Galaxy Catalogue - with a uniform isophotal detection limit 
of 26 mag arcsec -2 in the B-band - to revise the normalisa¬ 
tion parameters. This resulted in a reasonably consistent global 
Schechter function with: — 5 log h = (—19.64 ± 0.08) mag, 

<j>* = (1.64±0.005) x 10 -2 ft 3 Mpc -3 and a = -1.1±0.1 (taken 
from lPriver & De Pronrisf2003i) . 

Despit e the apparent convergence brought about by 
iLiske et ahl 1 2003 ) there are two strong reasons to now go 
beyond the monovariate Schechter function: 

(1) Numerical and semi-analytic simulations are sufficiently de¬ 
veloped that it is now possible to predict not just the luminos- 
ity di stribution but also the size distribution (see e.gjMc^yd] 
1 1 998t Buchalter, Jimenez & Kamionkowski 2001 : IShenet a "120031) 
and hence the bivariate distribution of luminosit y and size (or 
luminosity and surface brightness, see e.g. Sodre & Laha\||^93l 
Bo yce & Phillippslll995l : Ide Jong & Lacevl l200(1 ~ lcros!HiLDriven 
2Qo3h This bivariate brightness distri bution (BBD) can be fitted by 
a 6-parameter Choloniewski function 1 Choloniew skil 1 9851) - essen¬ 
tially a Schechter function with a luminosity dependent_Gaussian 
distribution in surface brightness (see e.g. Ide Jong & Lacevlhood : 
ICross & Drive J2 002) - or compared directly without confinement 
by any functional form. Either way, a multivariate distribution pro¬ 
vides a more stringe nt con straint than the monovariate Schechter 
fun ction (see al solDriveil2 004l) ; 

(2) ILiske et a DI3 essentially circumve nt the surface bright¬ 
ness issues raised by Cross & Drived 120021) by allowing the nor¬ 
malisation of the LFs to accommodate the missing flux from 
low surface brightness galaxies. This approximation succeeds 
because the galaxies that dominate the bright galaxy number 
counts are all around the L* value, where the surface brig htness 
distribution appears well defined (see iDrive 3 11999 . [Cross^eud] 

S . We also know from stud ies of giant ellipticals_^^^^^^ 

iGraham & Guzmadl1200 30, disk galaxies ide Jong <S^ace^ 


1 It is worth noting that the luminosity distributions seen in rich clus¬ 

ters such as A963, Coma and Virgo are also seen to tu rn up at a similar 

magnitude, (se e e.g. |Drivei^t^? ll 1994 iTrenthadl 199ijt Driver, C ouch & 
Phillipps 1998: lDrive^^^^ProDri^2QQ3r and iPopessoeta ll2004) t. Con- 
verse ly the local grou p lumino sity function remains flat to the Jeans mass 
limit iKarachentsev et all2003l) . 


2000), dwarf galaxies in the loc al group |Matedll993) , rich clus¬ 

ters iFereuson & Binggelll994t lAndreon&Culhan(Rd|200j)_and 
from early work on the global field BB D jprivej 1999 t]crossetaD 

l200lt IB lan ton et a Dbooit lsh en et a ilbooiT that there exists a lu¬ 

minosity surface brightness relation (see also I Bel 1 et aflEp OS for 
the near-IR relation and a discussion of its variation with wave¬ 
length). Together these imply that the impact of surface brightness 
selection effects will be a function of luminosity. Therefore recov¬ 
ering the full luminosity distribution can only be done by simulta¬ 
neously considering luminosity and surface brightness - a conclu- 
sion reached many times over the pa st half-century IZwic kyjl957|; 
|Pisnevll97fil : llmDev & Both uni 1 9971) . 

An additional motivation is that given the size of contempo¬ 
rary datasets, it becomes logical to start to explore multivariate dis¬ 
tributions to quantify trends such as the luminosity-surface bright- 
ness rela tion iDri veJ 19991) or the colour-lumino sity relation (e.g., 
iBaldrv et all2004l : see also lBlanton et all2003bl) . 

Here we reconstruct the bivariate luminosity-surface bright¬ 
ness distribution with careful attention to the impact of surface 
brightness and size selection effects. In Section^we describe the 
imaging and redshift data. In Section [3] we describe the conver¬ 
sion of our catalogue from apparent to absolute properties, the 
various selection boundaries and our methodology for deriving 
the joint luminosity-surface brightness distribution. We apply it 
to our data in Section |4] and compare our results to previous es - 
timates of the BBD JdTjo ng_&L i acevl|2QQ3; |Cross_etalJ |200l|), 
the LF iN orberg e t all2002l : 1zucca^Hln^jj^lBiantoi^Htll20(Hl 

iBlanton etall2 003a) and the surfac e brightness or size distribution 
iCross et alfcOOlUShen et alJ2003f) . We summarise our key results 
in Section[5] Throughout we assume H 0 = 0.3, Ha = 0.7 as de¬ 
fined in Section lTTI 


2 DATA 

2.1 Photometric catalogue 

The Millennium Galaxy Catalogue (MGC) is a deep (fj, UM = 
26 mag arcsec -2 ), wide-field (30.88 deg 2 ) B-band imaging survey 
obtained with the Wide Field Camera on the 2.5-m Isaac Newton 
Telescope. The survey region is a long, 35-arcmin wide strip along 
the equator, covering from 9 h 58 m 28 s to 14 h 46 m 45 s (J2000), and 
is fully contained within the regions of both the 2dFGRS and SDSS 
Data Releas e 1 ( DRhlAbazaiian et all2003}) . 

In lLiske et all 120031) we gave a detailed description of the ob¬ 
servations, reduction, object detection, Galactic extinction correc¬ 
tion and classification (using SExtractor: iBertin & Amout Jl99dh 
and presented precision galaxy number counts over the range 16 ^ 
B mgc < 24 mag. All non-stellar sources to B M oc = 20 mag were 
visually inspected and if an object was found to be incorrectly de- 
blended or if the object’s parameters were obviously wrong, the 
object was re-extracted by manually changing the SExtractor ex¬ 
traction parameters until a satisfactory result was achieved. In ad¬ 
dition, all low-quality regions in the survey (e.g. near CCD defects) 
were carefully masked out tsee ILiske et all2003l and lLemon et all 
2004). We demonstrated that the internal photometric accuracy 
of the MGC is ±0.023 mag and that the astrometric accuracy is 
±0.08 arcsec. In lCross et alJ t: 2004l) we then compared the photo¬ 
metric accuracy, completeness and contamination of the MGC to 
the 2dFGRS and SDSS-EDR and DR1 datasets. Being deeper and 
of higher resolution, we found that the MGC is both more accurate 
and more complete than either of these surveys. The MGC is cur¬ 
rently the highest quality and most complete representation of the 
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Figure 1. A comparison between Digital Sky Survey (left panels, bj-band), SDSS-DR1 (centre, g-band) and MGC (right panels) images: MGC03494 (upper); 
MGC18325 (middle); MGC04795 (lower). The data illustrate the greater depth and higher resolution of the MGC. All images are 33” on each side and scaled 
to comparable dynamic ranges.Figure degraded see http://www.eso.org/~jliske/mgc/ for full pdf copy of this paper. 


nearby galaxy population. Fig. Q shows images of three galaxies 
from the Digital Sky Survey (left), the SDSS (centre) and the MGC 
(right), illustrating the superior resolution and greater depth of the 
MGC data. 

The MGC-BRIGHT catalogue of iLiske et alJ hoQll ) was 
bounded at the bright end by 16 mag, where stars begin to flood 
in our best seeing fields. We now extend this catalogue bright- 
wards to B M ac = 13 mag by carefully examining all non-stellar 
objects in the range 13 ^ B M gc < 16 mag, repairing objects 
that had been over-deblended by SExtractor and checking com¬ 
pleteness with NED at B M gc < 15 mag, where SExtractor - op¬ 


timised for faint galaxy detection and photometry - may occasion¬ 
ally miss very large objects. A second change to MGC-BRIGHT 
with respect to iLiske et all |[ 2003l) is the deletion of some exclu¬ 
sion regions. Liske et al. had placed exclusion regions around very 
bright objects (_B M gc < 12.5 mag) to avoid spurious halo detec¬ 
tions and problems with background estimation in their vicinity. 20 
such regions had been erroneously placed around galaxies fainter 
than 12.5 mag and these have been deleted, resulting in an increase 
of the good-quality survey area from 30.84 to 30.88 deg 2 . 

Further minor changes were introduced to MGC-BRIGHT, 
each affecting ~tens of objects: different deblending decisions; 
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^best ^other 

Figure 2. Histogram of redshift differences for 2709 objects with mul¬ 
tiple redshift measurements. The distribution has an rms dispersion of 
92 km s — L The solid line is a Gaussian with the same median and rms 
as the data. 

spectroscopic identification of galaxies previously misclassified as 
stars (compact emission line galaxies) and vice versa (stars which 
appear extended because of almost perfect alignment with small, 
faint background galaxies); and identification of asteroids previ¬ 
ously misclassified as galaxies (discovered by comparison with 
SDSS-DR1 imaging data). 

The final MGC-BRIGHT sample now comprises 10 095 
galaxies. 

2.2 Redshifts 

2.2.1 Publicly available data 

By design the MGC survey region is fully contained within the 
2dFGRS and SDSS survey regions. Hence we first turned to these 
surveys in order to obtain redshifts for MGC-BRIGHT. In addi¬ 
tion we obtained redshifts from the 2QZ llCroom e~ alJ200lll . NED 
(excluding 2dFGRS and SDSS objects) and the smaller surveys 
of Francis, Nelson & Cutri (2004) (PF) and llmpev etal] 1 1 9961) 
(LSBG). In turn each survey, A, was matched against MGC galax¬ 
ies and stars with Bmoc < 22 mag. This magnitude limit and the 
following matching procedure were chosen in order to be able to 
account for every A' object and hence to verify the completeness of 
the MGC. 

First we found all pairs of A and MGC objects where each is 
identified as the other’s nearest neighbour in that catalogue. A pair 
with separation d smaller than some d cu t is classified as a match 
unless one or more other MGC objects have also identified the X 
object as their nearest neighbour in A with a separation < 2.5 d 
(where d refers to the separation of the closest pair). These X ob¬ 
jects as well as objects with d ^ d cu t and those which have not 
been identified as a nearest neighbour by any MGC objects are 
classified as initially unmatched. The value of d cu t is chosen as the 
separation where the invariably occurring tail in the distribution of 
separations begins to dominate. 

The above definition of initially unmatched objects is designed 
to conservatively flag all possible mismatches. Almost all mis¬ 
matches fall into one of three categories, which are a priori difficult 
to distinguish from one another but for which one needs to proceed 
differently: (i) unusually large positional disagreement; (ii) two (or 
more) objects in one catalogue where the other catalogue only has 


one object, because of different deblending choices or different res¬ 
olution; (iii) one object in X but none in MGC because of inaccu¬ 
rate photometry in either A or MGC, or because the X object is 
an asteroid or a spurious noise detection, or because the object is 
genuinely missing from the MGC. 

All initially unmatched objects are visually inspected in or¬ 
der to reliably identify those objects that were in fact correctly 
matched. The few cases where the mismatch was due to poor object 
extraction in the MGC were fixed. 

Each redshift of a matched object is assigned a redshift quality, 
Qz, which has the same meaning as in the 2dFGRS iColless et alJ 
1200II) : Q z = 1 means no redshift could be measured; Q z = 2 
means possible but doubtful redshift; Q z = 3 means probable red¬ 
shift; and Q z = 4 or 5 means good or very good redshift. Q z = 1 
and 2 are considered failures and we will use only redshifts with 
Qz ^ 3 in this paper. In addition to the 2dFGRS this system is 
also already used by PF. We translate the various quality param¬ 
eters of the other surveys to our system in the following manner: 
SDSS-DR1 redshifts are assigned Q z = 4 unless the SDSS ZSta- 
tus or ZWarning flags indicate a serious problem; for the 2QZ we 
use Q z = 5 —Q(2QZ); and NED and LSBG redshifts are assigned 
Qz = 3. 

In total the public data provided good quality redshifts for 
4804 MGC galaxies. 

2.2.2 MGCz 

In order to complete the redshifts for the full MGC-BRIGHT sam¬ 
ple we conducted our own redshift survey (which we label MGCz), 
targeting those galaxies which did not already have publicly avail¬ 
able redshifts. These observation were mainly carried out with the 
Two Degree Field (2dF) facility on the Anglo Australian Telescope. 
The data were collected, reduced and redshifted in an identical 
manner to that of the 2dFGRS (see iColless e~ all200ll for full de¬ 
tails). In particular we again used the same redshift quality param¬ 
eter, Q z , as the 2dFGRS (see above). 

The redshift incompleteness of the 2dF data showed a clear 
bias against low-surface brightness galaxies. Hence we under¬ 
took additional single-object long-slit observations, targeting low- 
surface brightness galaxies (/r® pp > 24 mag arcsec~ 2 , where fi^pp 
is the apparent effective surface brightness within the half-light ra¬ 
dius, see Section l3~4l as well as gaps in our 2dF coverage, with 
the Double Beam Spectrograph (DBS) on the Australian National 
University’s 2.3-m, the Low Resolution Spectrograph (LRS) on the 
3.6-m Telescopio Nazionale Galileo on La Palma, the ESO Multi- 
Mode Instrument (EMMI) on the New Technology Telescope in 
La Silla and the Gemini Multiple Object Spectrograph (GMOS, 
long-slit nod & shuffle mode) on Gemini North. In all cases we 
observed between 600 and 3600 s through a 1.5 to 2-arcsec wide 
slit, aligning the slit with the major axis of the object to collect 
as much flux as possible. All the data were reduced in a similar 
manner using standard IRAF procedures. The resulting spectra are 
of similar resolution as our 2dF data, but of higher quality: they 
increased our redshift completeness at B M ac < 19.5 mag and 
Mapp > 24 mag arcsec^ 2 from 90.7 to 98.8 per cent. 

2.2.3 Combined redshift sample 

Using all available data and MGCz data qthe overall redshift com¬ 
pleteness of MGC-BRIGHT is 96.05 per cent. In Fig.[3](third row 
from top) we show the final incompleteness as a function of appar¬ 
ent effective surface brightness and apparent magnitude, (g — r) and 
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Figure 3. Top row: the redshift incompleteness for the 10 095 galaxies in the 30.88 deg 2 MGC region using only 2dFGRS redshifts as a function of apparent 
effective surface brightness and apparent magnitude (left), (g — r ) (centre) and (u — r ) colour (right). The second row shows the same using only SDSS 
redshifts. These figures demonstrate the significant contributions of these surveys to our final redshift catalogue (third row) which has an overall completeness 
of 96.05 per cent. The light hashed region in the top left panel marks the approximate magnitude limit of the 2dFGRS, which varies as a function of position 
on the sky. The dashed lines in the second row mark the approximate SDSS spectroscopic magnitude and surface brightness limits. No galaxies beyond these 
respective limits were targeted by these surveys. The bottom row shows the number of objects per bin. 
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B „cc=13 mag 



Redshift (z) 

Figure 4. Absolute magnitude versus redshift for the MGC using the indi¬ 
vidual A'-corrections as defined in Section The solid curves show the 
bright and faint magnitude limits assuming the mean A’-correction and a 
universal evolution ofL = Ao(l-j-z) 0 ' 75 . The dashed lines again show the 
magnitude limits but now assuming the bluest and reddest A'-corrections 
used in the MGC. The upper and lower redshift limits are marked as vertical 
lines and the Mg value from Norberg et al. (2002) is shown as a horizontal 
dotted line. Note that the lower redshift limit defines the minimum absolute 
magnitude of Mg — 5 log h = —13 mag. The locations of known Abell 
clusters are indicated and clearly show up as redshift over-densities. 


(u — r) model colours. In TableQwe list the numbers of redshifts 
contributed by each of the surveys. 

For 2709 galaxies we have more than one Q z ^ 3 red¬ 
shift measurement, totaling 6226 redshifts. These duplications are 
mainly due to overlaps between the 2dFGRS and SDSS-DR1, but 
232 of these galaxies have at least one MGCz redshift. For these 
objects we pick a ‘best’ redshift by sorting first by Q z , then by 
signal-to-noise ratio (where available) and then by the order of the 
surveys in TableQ Fig.[2]shows the distribution of differences be¬ 
tween the multiple redshifts which has an rms of 92 km s -1 . 


3 COSMOLOGY, A -CORRECTIONS, EVOLUTION, 
SELECTION LIMITS AND ANALYSIS METHOD 

To calculate the MGC BBD we first need to adopt a cosmological 
framework, appropriate A'-corrections [A'(z)] and an evolutionary 
model [E{z)\. It is of course possible to determine a LF while 
simultaneously solving for one or more of these quantities - see 
e.g. lBlanton etaTJ i2Q03al) . who derive a LF while simultaneously 
solving for both the number-density evolution and the luminosity 
evolution. However, our concern with this approach is that it can 
lead to erroneous resu lts if the in comple teness is significant and 
biased (as indicated in Cross et ali l2004l and Fig. 0. In this case 
the evolutionary functions would be fitting not only evolution but 
also the uncorrected selection biases. Instead, our philosophy is to 
adopt reasonable assumptions for the evolution and concentrate on 
a strategy to account for selection bias. 




Figure 5. Two examples of our spectral fitting algorithm, indicating a good 
match to a young elliptical (upper panel) and older late-type (lower). Once 
the optimum spectrum has been determined this is used to determine the 
TAicc-band A'-correction. 

3.1 Cosmology 

We use Ho = 0.3 and Ha = 0.7 - broadly confirmed by the 
WMAP mission combined with the 2dFGRS power spectrum - but 
adopt h = fTo/(100 km s -1 Mpc -1 ) for ease of comparison with 
previous results. We also define maximum and minimum redshift 
limits. The maximum limit is chosen as the point at which red L*- 
galaxies fall below our imposed B M0C = 20 mag limit (this occurs 
at z M ax = 0.18, see Fig.0 excluding 1774 galaxies). The lower 
limit is determined by the point at which the redshifts are signif¬ 
icantly influenced by the local velocity field. From lPeacock et all 
<20011) we find that the pairwise velocity dispersion, as measured 
from the 2dFGRS, is 385 km s -1 . Hence to ensure minimal im¬ 
pact from peculiar velocities (i.e., Up ec „i ia r < O.ltWssion) we must 
set Zmin = 0.013 (excluding 44 galaxies). Fig.[4]shows the abso¬ 
lute magnitude versus redshift distribution for our data along with 
the upper and lower red shift limits_£vertical lines) and the typi¬ 
cal value of Mg iNorberg et alJl200a dotted line). Note that the 
lower z limit, combined with the faint apparent magnitude limit 
(-Bmcc = 20 mag) ultimately defines how far down the luminosity 
distribution one can probe (Mg — 5 log h « —13 mag). One can 
only reliably push fainter by extending the faint apparent magni¬ 
tude limit of the survey and not by covering a wider area (unless 
direct distances are measured or the detailed local peculiar veloc¬ 
ity field is modeled). In total there are 7878 galaxies within these 
redshift limits. 


3.2 The A-correction 

A-corrections, or bandpass corrections, are galaxy specific — 
ultimately galaxy component specific — and given the known 


variation in s 

Norberg et al. 

alaxy 

2002 

types, a mean A'-correction (as adopted in 
and in most previous studies) may be overly 

simplistic (see also 

Blanton et al. 2003a who reach a similar con- 
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Figure 6. The resulting ^^-probability distribution from our K(z) fitting 
algorithm (histogram) as compared to the expected distribution (for five 
degrees of freedom, solid line). Note that because we re-analyse galaxies 
with probabilities less than 5 per cent no values with xj? im ^11 occur. 

elusion). Here we identify a best fit synthetic spectrum from the 
combined broad-band MGC and SDSS-DR1 colours for each indi¬ 
vidual galaxy and then directly measure the A-correction from this 
optimal template. We use the extinction-corrected MGC Kron (2.5 
Kron radii) and extinction-corrected SDSS-DR1 Petrosian magni¬ 
tudes for this purpose. The Petrosian magnitudes are corrected to 
the AB system using u A b = Msdss — 0.04; g, r, i AB = g, r, i S Dss and 
~ab = 2sdss + 0.02 as recommended on the SDSS-DR2 website. 
The .Bmgc magnitudes are adjusted for the known zerop oint o ffset 
of 0.039 mag between the MGC and SDSS-DR1 l lCross et all2004l) 
and converted from Vega to AB. This conversion was found to be 
Bmgc(AB) = B M ac(Vega) — 0.118 for the MGC system. 

We now wish to derive the optimal spectral template match to 
these si x filte rs in the AB system. We elect to use the spectrum li¬ 
brary of lPoggiantil { 1997 ) which includes ellipticals and early and 
late-type spirals, with ages ranging from 15 Gyr to 3.7 Gyr (27 
synthetic templates in all). To determine the best template we calcu¬ 
late the AB magnitudes as if observed through the Bmgc and SDSS 
u,g,r,i,z filters for a single template and iteratively rescale the 
amplitude of the spectral template to minimise \ 2 • 

2 _ xp r Qc— smodel) i2 

/ '- £-^x=u,B ,g Air -I 

Where the errors are given by: Au = ±0.10; AB = ±0.07; A g = 
±0.05; Ar = ±0.05; A i = ±0.05 and Az = ±0.08 mag. These 
error values were derived as follows. Firstly, using just the SDSS 
filters, we determined the xiLrdistribution for the optimal spec¬ 
tral fits based on the errors listed on the SDSS website (Au = 
±0.1 [red leak], A gri = ±0.02, A z = ±0.03). We then scaled 
the griz errors until the recovered ^^-distribution matched the ex¬ 
pected distribution for four degrees of freedom. We then repeated 
this process now including the B M ac filter (see Fig.|6j. These final 
errors are in good agreement with the study of lCrosset afl hoQdt) . 

The \ 2 -minimisation was repeated for each of the 27 spec¬ 
tral templates and the smallest of the 27 xiL values used to iden¬ 



Redshift (z) 

Figure 7. Our final /7-corrections. The MGC 17-corrections trace out the 27 
tracks of the synthetic spectra of Poggianti (1997). The filled circles indicate 
the mean of our data and the vertical dashed lines show the adopted redshift 
limits used in the BBD analysis. 


tify the most appropriate template. Finally if the xiL "Probability 
for the optimal template was less than 5 per cent we repeated the 
process, first rejecting the S MOC filter and then the u and z filters 
in turn (hence the truncation at xiL ~ 11 on Fig.|6j. Magnitudes 
fainter than the quoted 95 per cent completeness limit were not 
use d (u = 22.0, q = 22 .2, r = 22.2, i = 21.3, z = 20.5 mag, 
see lAbazaiian et all2003t) . Fig .|7|shows two example galaxies with 
the best-fit spectral templates overlaid on the measured magnitudes. 
Fig .[6] shows the resulting xlL distribution (histogram) and the ex¬ 
pected x' 2 -distribution (for five degrees of freedom, solid line). 

Shown in Fig.0are the resulting A'-corrections for the MGC, 
which trace out the 27 different spectral templates, as well as the 
mean A'-correction in 0.01 redshift intervals (large dots). Finally, 
in the few cases (25 galaxies in total) where a galaxy did not have 
an SDSS-DR1 match we assigned the mean A-correction value at 
the galaxy’s redshift. 


3.3 Evolution 

We adopt a global evolution of the form L = Ao(l ± z) 13 with an 
ini tial value of 0 = 0.75 (i .e., E(z) = —0.75 x 2.51og(l ± z), 
see tPhillipDS & Driver! 199.‘f) . Given the redshift limits this implies 
an evolutionary correction in the range 0.01 — 0.14 mag. This is a 
passive luminosity evolution model with no merging over this red¬ 
shift range (0.013 < z < 0.18). In Section l4~3l we solve for /? and 
find —2.0 < j3 < 1.25 (see Fig. 1201. More complex evolution¬ 
ary scenarioes will be considered in future papers where we sub¬ 
divide the galaxy population by t ype and com ponent. However we 
do n ote that recent studies (e.g.. lPatton et alJl200(il) including our 
own toe Propris et alJ2005h indicate a very low low-z merger rate 
based on counts of dynamically close pairs. It is also worth noting 
that the step-wise maximum likelihood method used assumes no 
number density evolution. 
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3.4 Estimating MGC effective surface brightnesses 

We adopt the effective surface brightness, i.e. the average surface 
brightness within the half-light radius, as our structural measure¬ 
ment. The effective surface brightness is chosen as it makes no as¬ 
sumption of the radial flux profile and can be measured empirically 
from the data for any galaxy shape. Hence it is a direct measure 
reflecting the global compactness regardless of whether the galaxy 
is smooth or lumpy. 

The half-light radius is the semi-major axis of the ellipse con¬ 
taining half the flux of SExtractor’s BEST magnitude. The ellipse’s 
centre, ellipticity and position angle are taken from SExtractor. 
For compact objects the half-light radius is affected by the see¬ 
ing. Through simulations of both exponential and de Vaucouleurs 
profile galaxies using IRAF ARTDATA (see Appendix[^]for details) 
we find that we can correct the observed half-light radius for this 
blurring and recover its true value, r G , to an accuracy of 10 per cent 
(see Fig. lA2t by: 

r° = vV|-0.32r 2 , (1) 

where r e is the observed half-light radius and T is the FWHM of 
the seeing profile, all in arcsec. 

Work is currently underway to determine 2D bulge-disk de¬ 
compositions from the MGC data which will be presented in a fu¬ 
ture paper. However at this stage it is unclear whether the effective 
radius or independent bulge and disk parameters are the more fun¬ 
damental. Here, we define our absolute effective surface brightness, 
fi e , as: 

M e = Hmgc + 2.5 log[27r(r“) 2 ] -10 log(l + z) - K(z) -E(z), (2) 

We note that r G is defined along the semi-major axis and hence the 
above prescription assumes galaxies are optically thin, i.e., derived 
parameters are independent of inclination. The role of dust will be 
explored in a future paper. 


3.5 Selection limits 



Figure 9. The distribution of galaxies with (dots) and without (open trian¬ 
gles) redshifts along with selection limits. Square data points with errorbars 
show the selection limits as determined from the simulations of Appendix^] 
and the solid line shows the fit. The low surface brightness and maximum 
size limits are defined as limits of reliable photometry by requiring errors in 
the photometry of less than ±0.1 mag and ±0.1 mag arcsec -2 . The dashed 
line is the fundamental detection limit. 

fective surface brightness plane, along with the galaxy distribution 
with (dots) and without (triangles) redshifts. 

For each galaxy i with known redshift, Zi , we can now deter¬ 
mine the parameter space in absolute magnitude, M , and absolute 
effective surface brightness, fi e , over which this galaxy could have 
been observed, given the survey limits. The appropriate limits are 
(following lDriveJ1999t) : 


Like all surveys our catalogue will suffer from some level of in¬ 
completeness in both the imaging and spectroscopic surveys. As 
our data are the deepest available over this region of sky we cannot 
empir ically quantify our imaging incompleteness (see ICross et all 
12004 for an assessment of the incompleteness of the SuperCOS- 
MOS Sky Survey, 2dFGRS and SDSS-EDR/DR1 relative to the 
MGC). We can, however, define the selection limits at any red- 
shift, following lDrive r 11999), and use this information later (Sec- 
tion l3.6> to quantify the selection boundary in our final space den¬ 
sity distribution. 

We begin by noting the MGC’s apparent survey limits: a 
bright and faint apparent magnitude limit of m BRIGHT = 13 mag 
and m FA i NT = 20 mag, a limiting apparent effective surface bright¬ 
ness of /r RIM = 25.25 mag arcsec -2 , a minimum size of r Mm = 
(0.04T 2 + 0.37T + 0.1)5 arcsec (where T is the seeing FWHM) 
and a maximum size of r MAX = 15 arcsec. The last three limits 
refer to the seeing-corrected quantities (see Section l3~4l and were 
derived from simulations of 45°-inclined optically thin exponen¬ 
tial disks using IRAF’s ARTDATA package (see Appendix[X]for de¬ 
tails). It is important to note that the surface brightness and maxi¬ 
mum size limits do not define the absolute detection limits but the 
boundary of reliable photometry (where ‘reliable’ is arbitrarily de¬ 
fined as ±0.1 mag and ±0.1 mag arcsec -2 accuracy, see Fig. lA2> . 
Fig. |9| shows these limits in the apparent magnitude-apparent ef¬ 


M 


BRIGHT,! 

TT faint? i 

e 

/h ow, i 

e 

Mlow,w 


= m BRIGHT - 5 log d L (zi) - 25 - K(zi) - E(zi) 

= m FAINT - 5log d L (zi) - 25 - I<(zi) - E(zi) 

= Mum - 101og(l + Zi) - K(Zi) - E(Zi) 

= Mi+ 5 log d L ( Zi) + 25 + 2.5 log(27rr 2 AX ) (3) 

— 101og(l + Zi) 

= Mi+ 5 log d L (Zi) + 25 + 2.5 log(27rr 2 N ) 

— 101og(l + Zi) 


Note that )z BOW t is defined by both the minimum effective surface 
brightness for reliable detection, /r BIM , and the maximum reliable 
size limit, r MAX . The maximum reliable size limit is imposed by the 
smoothing box size when calculating the bac kground sky levels. 
These limits are introduced in more detail in iDriverl dl999 ). One 
can also envisage a sixth limit (also /z® G „ ;i ) due to the dynamic 
range of the CCD detector. In our case none of the galaxies have 
their central regions flooded and so this limit can be ignored. 

Fig. llOl shows the selection limits for three galaxies at redshifts 
0.025 (lower panel), 0.100 (middle) and 0.175 (upper) along with 
the location of the galaxy (large solid circle), galaxies with similar 
redshifts (open squares) and the full (0.013 < z < 0.18) MGC 
distribution (dots). The selection boundary consistently forms a 
five-sided figure which glides through the full distribution as the 
redshift increases. Because of the individual iC-corrections each 
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Figure 8. An illustration of our BBD SWM1 approach. Left panels: examples of the individual matrices U,, (weighting function; upper) and /lijk (pseudo¬ 
visibility function; lower) for MGC16504. Right panels: Wijk (upper) and ^ IIjj (lower).Figure degraded see http://www.eso.org/~jliske/mgc/ for 

full pdf copy of this paper. 


galaxy has its own unique selection window. It is worth noting 
that large luminous galaxies are undetectable, or at least photo¬ 
metrically unreliable, at very low redshift. Hence the methods of 
IShen et al . 1 2003 ) which use z MEmAN (M, /r e )-style corrections, may 
underestimate this population because they implicitly assume se¬ 
lection boundaries only cut into the distribution at some upper red- 
shift limit. lCross et al J b 2()0 1 ) accounted for this effect in their anal¬ 
ysis. 

For each galaxy with known redshift we now define an ob¬ 
servable window function, Oi{M, /r e ), as: 

{ 1 if M B8IGHT ,i < M < M FAmT ,i 

and Uman,i < M e < Plow,* (4 ) 

0 otherwise. 

The functions Oi will be used in Section lT6l in the construction of 
the BBD. 


3.6 Bivariate step-wise maximum likelihood 

We now construct the MGC BBD using a bivariate brightness step¬ 
wise maximum likelihood (SWML) estimator similar to that de¬ 
scribed bv ISodre & Lahav 1 ( 11993 ). This is an extension of the stan¬ 


dard SWML method defined bv lEfstathiou et al .. ( 1988 ) in the sense 
that the sample is now divided into bins of both absolute magni¬ 
tude and absolute surface brigh tness as opposed to just absolute 
magnitude alone. Note that lSodre & Lahav . ( 1993 ) formulated their 
method in terms of absolute magnitude and a physical diameter. 
However, the MGC is primarily magnitude and surface brightness 
limited with a fixed isophotal detection limit of 26 mag arcsec^ 2 . 
For this reason we elect to work in M-jT space rather than M-R e 
space. 

Essentially, for i = 1...N objects the volume-corrected rel¬ 
ative BBD in j = 1... J absolute magnitude bins and k = 1...K 
surface brig htness b ins with widths AM and A/j, e can be evaluated 
by (cf. Sodre & Lahavl 19931) : 


— 


Ef Wyfc 


Ef[#W Ef E™ 4>lmHilm] ’ 


(5) 


where the relative space density of galaxies, <f>jk, is given in terms 
of the space density of the previous iteration, the weighting 
function, Wijk, and the visibility function, Hijk- 

The weighting function, Wijk , corrects for the incomplete¬ 
ness in our redshift survey. This incompleteness is demonstrably 
not random (see Fig.0 and we are preferentially missing faint dim 
galaxies with extreme colours. When constructing LFs this bias is 
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MGC27516 at 2=0.175 (0.174 < z < 0.176) 



-22 -21 -20 -19 -18 -17 -16 -15 -14 


Absolute magnitude (M B , mag) 
MGC17731 at z=0.1000 (0.099 < z < 0.101) 




Absolute magnitude (M B , mag) 


Figure 10. The selection boundaries for three galaxies as indicated along 
with the location of the galaxy (solid circle), those with similar redshifts 
(open squares but differing if-corrections) and the full (0.013 < z < 
0.18) MGC distribution (dots). The figure illustrates that the selection 
boundaries vary dramatically with redshift (bottom to top) such that each 
redshift slice, within the survey limitations, only samples a restricted region 
of the luminosity-surface brightness plane. Note that some points may lie 
outside the box if their If-corrections are significantly different from the 
galaxy for which the selection lines have been drawn. Figure degraded see 
http://www.eso.org/~jliske/mgc/ for full pdf copy of this paper. 



Figure 11. The normalised bivariate brightness distribution. The vertical 
axis is logarithmic in units of h 3 Mpc -3 mag -1 (mag arcsec -2 ) -1 . 


Davies & Disney 1990), although not precisely as it does not ac¬ 
count for the variation of physical cross-section with redshift or 
for isophotal corrections. The steep sides of the Hijk surface 
(Fig- ID highlights the narrow range over which surface brightness 
selection effects go from negligible to extreme and hence the im¬ 
portance of understanding a survey’s selection limits. 


often ignored (e.g.. lBlanton_et_al] |2Q03a|) or considered a function 
of magnitude only (e.g.. lNorbe^e^il20Q2l) . Here, we account for 
luminosity and surface brightness dependent incompleteness and 
define Wijk as: 




if M ■ — AM < A/f. <- M _l_ AM 
Ni(Q z ^ 3) 11 1V1 0 2 ^ iV1% ^ 1V1 3 ^ 2 

and nl - < nl < ,j,l + *£. (6) 

0 otherwise, 


where Ni is the total number of galaxies lying in the same ap¬ 
parent magnitude-apparent surface brightness bin as galaxy i and 
Ni(Q z ^ 3) is the number of galaxies with known redshifts (i.e., 
those with Q z ^ 3, see Section |3 in the same bin. Hence each 
galaxy is weighted by the inverse of the redshift completeness in 
the galaxy’s apparent magnitude-surface brightness bin. 

The visibility function Hijk is given by: 


(AMA(i e ) 


/ AM 

i Mj —AM/2 


d M e Oi(M, M e ),(7) 

uf-A^/2 


where Oi(M, ji e ) is defined in Section l3~5l Hence Hijk is the frac¬ 
tion of the Mj-fj,\ bin which lies inside the observable window of 
galaxy i. 

Fig. 0 shows example Wijk and Hijk matrices for i = 
MGC16504 (upper and lower left respectively) as well as the 
summed IV,,,y, and Hijk matrices (upper right and 

lower right respectively). Wijk is essentially the observed 
galaxy distribution corrected for incompleteness and ^. Hijk is 
analogous to the visibility surface of Visibility Theory (Phillipps, 


4 THE MGC BIVARIATE BRIGHTNESS DISTRIBUTION 

We now apply the above methodology to our MGC sample. To re¬ 
cap, we extract galaxies in the range 13 < B Mac < 20 mag, de¬ 
termine individual A'-corrections, determine seeing-corrected sur¬ 
face brightness estimates, attach an incompleteness weight to each 
galaxy, reject galaxies which lie outside of their designated reli¬ 
ably observable regions (7 galaxies too large, 99 too compact and 
22 too dim), generate weighting and visibility matrices for each 
galaxy and then iteratively apply equation J5J until a stable set of 
(p(M,[i e )-v alues is found. Fig. IIII shows the result of this proce¬ 
dure after renormalisation (described below). It is worth noting that 
the majority of the rejected galaxies are extremely compact and 
have entered our sample only through spectroscopy of apparently 
stellar objects (2QZ, SDSS QSO survey, NED and MGCz). 

The SWML method, by its construction, loses all information 
regarding the absolute normalisation (cf. lEfstathiou et alJft98 8). 
Here the normalisation issues are compounded by the considera¬ 
tion of selection bias. We overcome this by calculating the num¬ 
ber of galaxies in the absolute magnitude and absolute surface 
brightness interval —20 < Mb — 51ogh < —19.5 mag and 
20 < fj, e < 24 mag arcsec -2 . This parameter space is observ¬ 
able for any K(z) over the redshift range, 0.047 < z < 0.162 
(equivalent to a co-moving volume of 311849 b -3 Mpc 3 for the 
30.883 deg 2 field-of-view). Within these redshift, absolute mag¬ 
nitude and surface brightness ranges there are 785 galaxies with 
a summed weight, due to incompleteness, of 787.8. The <j> val¬ 
ues are now rescaled to reproduce the correct space density for 
these galaxies (i.e., 0.00253 ± 0.00020 h 3 Mpc -3 ). The error in- 
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eludes both Possion and large scale structure considerations and 
was derived from mock 2dF NGP catalogues t Cole et al 1998 : 
http://star-www.dur.ac.uk/~cole/mocks/main.html i extracted from 
the Virgo Consortium Hubble Volume. Volumes equivalent in shape 
to the MGC normalisation volume were extracted and the standard 
deviation of —20 < Mb — 5 log ft < —19.5 mag galaxies deter¬ 
mined. The final error, inclusive of large scale structure, is roughly 
twice that of the Poisson error alone. 

Fig. 1121 shows the final normalised BBD on a logarithmic 
scale. The thick dotted line shows the region within which the 
statistical errors are 25 per cent or less. The thick solid line 
shows the effective detection limit defined as the BBD region sam¬ 
pled by at least 100 galaxies (equivalent to a volume limit of 
~ 1800 ft -3 Mpc 3 ). The detection limit bounds the BBD at the 
low surface brightness and faint limits and clearly starts to im¬ 
pinge on the distribution from around Mb — 5 log ft « —18 mag 
(coincidentally the usually adopted boundary between the dwarf 
and giant galaxy populations). This confirms that even the MGC, 
the deepest local imaging survey to date, remains incomplete 
for both high {Mb — 5 log ft sa —15 mag) and low surface 
brightness galaxies. To sample these regions requires both deeper 
(Mlim 26 mag arcsec -2 ) and higher resolution imaging data 
(T 1.27 arcsec) over a comparable or larger area (Q 2? 30 deg 2 ) 
with high spectr oscopic completen ess. 

Following Jone & Lacexl fcood ) and [c ross & Privet! 
120021) we attempted to fit th e MGC BBD with a Choloniewski 
function lChotoniewskill985 ). given by: 


= °ll nl 0 <j>* 10 O - 4 ( M * -M)(o,+l) e _ to°- 4 ( M * ~ M ) 

V27r oy,e 


Xexp -- 




( 8 ) 


The parameters < />*, M * and a are the conventional Schechter func¬ 
tion parameters J.Schechteil 197fil ). The three additional parameters, 
p e *, and /3 are the characteristic surface brightness, the Gaus¬ 
sian width of the surface brightness distribution and the luminosity- 
surface brightness relation respectively. This f unction has been 
deri ved analytically b y iDal canton et all j 19971) . iMo et alJ 1 1 9981) 
an d|de Jong & Lacev (2000). Note that the additional normalisa¬ 
tion factor, (v / 27rcr f i<0 _1 > ensures that <f >*, M* and a are directly 
equivalent to the more familiar Schechter function parameters. The 
fit is achieved using the downward simplex method ((Press et al1 
11992) and the 1 -a errors are derived from a one-dimensional y 2 - 
minimisation across each parameter (ltr = Ay 2 = 158.9). 

For the 151 degrees of freedom (157 data values minus 6 
model parameters) we find an unreduced y 2 ,, of 608.7. This sug¬ 
gests that the Choloniewksi function is a poor fit to the data (see 
Fig. 1 1 3t . The main reason appears to be the intrinsic assumption of 
a constant cr M e, whereas the distribution show n in Fig. IIII clearly 
broadens, a feature also observed by n en et alJ j2003i) in their 
SDSS data. We discuss this further in Section l4~2l Choloniewski 
parameters and associated errors are tabulated in the upper half 
of Tabled along with other published data inclu ding outyprevious 
measurement based on the 2dFGRS iCross & Driver«2002t) . 


4.1 The space density of galaxies 

Of particular interest over the past few decades has been the galaxy 
luminosity distribution. This can be derived from the BBD by inte¬ 
grating over the surface brightness axis, (the associated errors are 
the root sum squares of the individual errors for each bin). Fig. 1141 


-22 -20 -18 -16 -14 

MB UGC -5 |0 9 h| (mag) 



Figure 12. The final space density or bivariate brightness distribution of 
galaxies shown as both a greyscale image and contours on logarithmic 
scales in units of ft 3 Mpc -3 mag -1 (mag arcsec -2 ) -1 . The contour spac¬ 
ing is 1 dex in O with the greyscale at 0.5 dex intervals. The thick solid line 
denotes the selection boundary defined as the BBD region where at least 
100 galaxies could have been detected and is equivalent to an effective vol¬ 
ume limit of 1800ft -3 Mpc 3 . The thick dashed line encompasses the region 
within which the statistical error is smaller than 25 per cent. 


(left panel) shows the result: <j>* = (0.0177 ± 0.0015) ft 3 Mpc 3 , 
Mg — 5 log ft = (-19.60 ± 0.04) mag and a = -1.13 ± 0.02. 
Note that the error in the normalisation parameter, <J>*, is the combi¬ 
nation (in quadrature) of the error from the fitting algorithm com¬ 
bined with the error from the number of galaxies within the nor¬ 
malisation volume. For completeness we also show the implied LF 
from our Choloniewski function fit (dotted line). However recall the 
poor Xmin value from this fit which highlighted the fit’s inability to 
model the dwarf population correctly (mainly due to the broaden¬ 
ing of the surface brightness distribution). The right panel of Fig fT-il 
shows the formal 1,2,3-cr error contours for the Schechter function 
fit to the collapsed MGC BBD luminosity distribution. 

4.1.1 Comparison to the ESP, 2dFGRS and SDSS LFs 

Fig. ll5l compares our LF based on the BBD analysis of the deeper, 
higher quality, higher completeness MGC data to a number of other 
recent LF estimates (tabulated in the lower half of Tabled- Note 
that these Schechter functions are sh own with their original cf>*- 
values as opposed to the iLiske et alJ 1200 3) corrected values. We 
find our LF recovers a comparable M* value to most previous 
surveys, but a slightly higher normalisation (as one might expect 
from a BBD-style analysis which recovers more light from the low 
surface brightness population than a non-BBD style analysis), and 
an a value consistent with the median of these other surveys. In 
comparison to the ESP LF iZucca et all 19971) we find consistent 
Schechter function parameters within the quoted errors of the two 
surveys. In comparison to the 2dFGRS LF JNorberg et all2Qo3) we 
find consistent M* and <p * values but a flatter (move positive) a. 
This discrepancy in a is formally significant at the 3cr-level, how¬ 
ever we do note that the 2dFGRS has undergone a major overhaul 
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Figure 13. The contours show the Chotoniewski function which best fits our 
final BBD in units of /i 3 Mpc —3 mag -1 (mag arcsec -2 ) —1 . The greyscale 
image and the thick solid line are the same as those in Fig. El 


Figure 15. A comparison of recent luminosity function estimates as indi¬ 
cated (left panel) and the associated errors in M* and a together with the 
MGC 1, 2 and 3-cr contours (right panel). 



Figure 14. Left panel: the luminosity function derived by integrating 
the BBD over surface brightness with the best-fit Schechter function 
(solid line). Also shown is the implied Schechter function from the fitted 
Chotoniewski function (dashed line). Right panel: the 1, 2 and 3-cr contours 
of the Schechter function fit. 


o f its photometry since this result (see discussion in ICross et al J 
120041) . The final 2dFGRS post-recalibration LF has yet to be de¬ 
termined. Our result are consistent with the M* value only for 
SDSS1 and inconsistent with all three Schechter function values 
for SDSS2. A detailed discussion of this discrepancy follows in the 
next section. 


4.1.2 Exploring the discrepancy with the SDSS LFs 

The most notic eable outliers in Fi g. 1151 are the t wo SDSS re- 
suits SPSS 1 telanton et alJ 1200 ll) and SDSS2 JBlanton et alJ 
12003 al) . The discrepancy in the nor malisation of SDSS1 has 
been comprehensively discussed by both N orberg et alJl2002l) and 
iBlanton et alJ 12003 al) and ascribed to the lack of an evolutionary 
model in the SDSS 1 analysis, coupled with their method of normal¬ 
isation. An over-density in the SDSS commissioning data region 
alsc^contributes. The SDSS1 normalisation has been revised by 
iLiske et °aT] 120031) : </> s *dssi(NEW) = 0.76 (OLD). Once cor¬ 
rected, the SDSS1 LF agrees well with the MGC (2<r-discrepancy 
in a), the 2dFGRS and the ESO Slice Project (ESP) estimates. 

The more recent SDSS2 result has a significantly fainter M * 2 
and a much flatter faint-end slope than any of the other surveys. 
|B lan ton et alJ l2003al) argued that the difference in the r-band LFs 
of SDSS1 and SDSS2 was simply due to the inclusion of luminos¬ 
ity evolution in the latter and that the luminosity densities measured 
by the 2dFGRS and SDSS2 agreed to within the errors once the dif¬ 
fering amounts of evolution had been taken into account. However, 
in the g-band the discrepancy in a between SDSS1 and SDSS2 
is more than twice that in the r-band. Furthermore, even though 
the luminosity evolution derived in the SDSS2 fitting process is 
very strong (Q = 2.04 mag per unit redshift, where M*(z) = 
M* (z = 0) — Q z) it cannot explain away the difference between 
the MGC and SDSS2 values. Our evolutionary correction given in 
Section l3~3l corresponds to Q « 0.74. If we replace this with the 
SDSS2 evolution, E(z ) = —2.04 a, and re-derive our collapsed 
MGC LF we find Mg — 5 log ft = —19.36 mag and a = —1.08 
(or Mg — 5 log ft = —19.36 mag and a = —1.06 if we restrict the 
LF fitting to the equivalent SDSS2 limit of Mg — 5 log ft < —16 
mag). Hence some tension remains between the SDSS2 and MGC 

2 Note that the conversion from 01 g to Bmgc is given bv lBlanton et all 
•2003a) as S M cc = 01 g - 0.09. 
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Figure 16. The left panel shows the result of a re-analysis of the MGC data 
with various low surface brightness selection limits imposed (as indicated). 
This figure demonstrates the extreme sensitivity of the recovered Schechter 
function parameters to the spectroscopic surface brightness completeness 
limit. The right panel shows 1 -a error ellipses for the Schechter function 
fits along with the data points from the right panel of Fig.L15J 


LFs irrespective of evolution. The same argument holds when com¬ 
paring SDSS2 to the 2dFGRS since the latter adopted an evolution 
equivalent to Q sa 1.15. 

To test whether the surface brightness selection limits might 
be responsible. Fig. 1 161 shows a re-analysis of the MGC data, fol¬ 
lowing the same procedures as before, but now imposing progres¬ 
sively shallower surface brightness limits (as indicated). The left 
panel again shows the LFs and the right panel the 1-er error ellipses 
of the Schechter function fits. It is quite apparent that shallower 
surveys progressively miss a greater fraction of the dwarf popula- 
tion(s) - and ultimately giants - recovering both a fainter M* and 
a more positive a. In fact, Fig. ll6l based on MGC d ata m irrors the 
expectation from the simulations of lCross & Driven ( 20021 ) remark¬ 
ably well (see their Fig. 5). 

On Fig. the LF with an imposed selection limit of 
23.25 mag arcsec -2 appears to match the SDSS2 LF closely and 
initially suggests that surface brightness selection might explain 
the SDSS2 result. (Note that SDSS1 implemented a BBD-style 
analysis). However Fig. 0 (left panels) do not show a signifi¬ 
cant difference in the low surface brightness spectroscopic com¬ 
pleteness between the 2dFGRS and SDSS surveys suggesting lim¬ 
its closer to 24.25 mag arcsec -2 . We note that the SDSS2 sur¬ 
face brightness 50 per cent completeness limit (including con¬ 
sideration of deblending and flux loss is sues , etc.) is_c|uoted as 
/r®[Circular] « 23.4 mag arcsec -2 felanton et alJl20053 . their 
fig. 3), where [Circular] indicates that the SDSS surface brightness 
measures are derived from ‘circular’ 3 half-light radii. Conversely 

3 The SDSS2 surface brightness measure is derived from the combination 
of SDSS Petrosian magnitudes and Sersic half-light radii based on circular 
aperture photometry. At this point it is not clear why the SDSS team have 
adopted circular apertures and how one should interpret circular half-light 


the MGC (and most contemporary datasets) use major axis half- 
light radii. These are derived by growing the optimal ellipse until it 
encloses 50 per cent of the flux. The use of a circular aperture will 
always result in a lower half-light radius than an elliptical aperture 
(except for purely face-on systems), and hence in a higher effec¬ 
tive surface brightness. The two effective surface brightnesses are 
related by: 


/T [Major Axis] 


/i® [Circular] — 2.5 log 



(9) 


where ( b/a ) is the axis ratio. For a random distribution of infinitely 
thin disks the expected median axis ratio (b/a) is 0.5. In real¬ 
ity the galaxy distribution is more diverse consisting of spheroids, 
bulges and thick disk components/systems. Circular half-light radii 
measurements are also somewhat resolution dependent, introduc¬ 
ing both seeing and redshift dependencies into the conversion. Em¬ 
pirically we find for our MGC data: (b/a) mgc = 0.71, imply¬ 
ing a correction of 0.37 mag arcsec -2 . Although we cannot di¬ 
rectly compare our MGC elliptical half-light radii to the circular 
Sersic half-light radii used by SDSS2, we can compare them to 
the standard SDSS-DR1 circular Petrosian half-light radii. We find 
(t*50,Petrosian/fMGc) = 0.84, independently confirming the cor¬ 
rection of 0.37 mag arcsec -2 . Hence the quoted SDSS2 50 per cent 
completeness limit of )t® [Circular] = 23.4 mag arcsec -2 equates 
to /r® [Major Axis] = 23.77 mag arcsec -2 . 

The colour correction is more complex as it depends strongly 
on galaxy type and luminosity. The low luminosity population is 
the bluest ((B M ac — r) ~ 0.5 mag) and hence surface bright¬ 
ness incompleteness will preferentially affect the faint end of the 
LF. Evidence for such a colour completeness bias can be seen in 
the observed SDSS incompleteness within the MGC region. This 
is shown as a function of (g — r) or (u — r) versus ulpp in 
the middle row (centre and right panels) of Fig. 0 As the SDSS 
spectroscopic survey is r-band selected it should not be surpris¬ 
ing that the B M oc completeness is a more complex surface in mag¬ 
nitude, surface brightness and colour. In the worst case scenario 
we find a 50 per cent completeness limit for the bluest systems of 
Ub MC iC = 24.27 mag arcsec -2 which is consistent with Fig.0( mid¬ 
dle left), but more importantly significantly fainter than the surface 
brightness limit implied by Fig.^| That a colour selection is obvi¬ 
ous in both the ( g — r) and (u — r ) SDSS completeness maps of 
Fig.^argues that the colour selection bias is not a feature of the u 
band data but more a general trend across all filters. 

Adopting an appropriate Ms mgc selection criterion is clearly 
non-trivial, however a simple way forward is to construct a 
LF from the available SDSS2 redshifts within the MGC re¬ 
gion (see Table. 1). This will have the SDSS selection func¬ 
tion inherently built in. We adopt m FAINT = 18.5 mag, /r® IM = 
24.27 mag arcsec -2 , the SDSS2 g -band evolutionary correction 
and apply a conv entional monovariate step-wise maximum likeli¬ 
hood fefstathiou et alll988l) . We now recover Mg = (—19.26 ± 
0.07) mag and a = -0.83±0.06 (or Mg = (-19.26i0.07) mag 
and a = —0.82 ± 0.06 for Mg < —16 mag) and consistent at the 
~ 1-cr level with SDSS2. We therefore conclude that the discrep¬ 
ancy between the MGC and SDSS2 LFs is complex, but predom- 


radii for moderate and highly inclined systems, particularly given the ad¬ 
ditional dependencies this will introduc e as a function of both seeing and 
distance. It is also not clear why|Blanton_etaU|2003a^_grefer_Sersic fitted 
half-light radii rather than their more reliahkT mlantoTct all2005bh empiri¬ 
cally measured Petrosian radii (which agree closely with MGC circularised 
half-light radii). 
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inantly due to the colour bias within the SDSS (cf. Fig. 0 middle 
right). This highlights the additional complexity in measuring LFs 
for filters other than the spectroscopic selection filter. 


4.1.3 Lower limits on the very faint-end 

Finally we note the recent preprint by iBlanton et all fcOOfiah in 
which an attempt is made to sample the very low luminosity 
regime from their New York University Value Added Catalogue 
IBlanton et alJl2005bl) . This incorporates a strategy for correcting 
the surface brightness incompleteness within the SDSS2 dataset 
resulting in a significantly steeper faint-end slope (a = —1.28 to 
— 1.52 in r). The premise is that the surface brightness distribution 
of the dwarf population can be extrapolated from that of the gi¬ 
ant population. This results in correction f actors for the spa ce den¬ 
sity of dwarf systems of up to 3.5 (see IBlanton et alJl2005al their 
fig. 6). In a similar spirit of producing a ‘best guess’ LF at low 
luminosities Fig. ll7l shows the MGC LF with the selection bound¬ 
aries and lower redshift bounds removed. This undoubtedly repre¬ 
sents an underestimate of the true space density at low luminosities 
(Mb — 5 log h > —14 mag), but in so far a s a c ompa rison can be 
made it is consistent with the SDSS result IBlanton et a l]2005al) . in 
the sense that the SDSS data lie above the MGC lower limits. Our 
LF is also consistent with a compilation of 231 galaxies within the 
local sphere (D < 5 Mpc) taken from I Karachentsev ~ all <2004 . 
Apart from indicating the amount of work that is yet to be done, 
Fig. ll7l does suggest that the low redshift global luminosity func¬ 
tion must steepen beyond a = —1.13 at some intermediate magni¬ 
tude (Mb — 5 log h > —16 mag). However it would seem unlikely 
that it could steepen enough for the dwarf systems to make a signif¬ 
icant contribution to the total local luminosity density (as indicated 
by the dashed line). 



Figure 17. Tentative constraints on the very faint end of the galaxy luminos¬ 
ity function from three independent studies. Solid squares: available MGC 
data using the BBD-style analysis but with all limits removed. These data 
are incomplete at the faint end with the loss of both high and low surface 
brightness systems and therefore repres ent lower limits only . Open circles: 
a recent ‘speculative’ SDSS estimate IBlanton et alJl2005al) based on an 
extrapolation of the surface brightness distribution at bright magnitudes. 
Crosses: Local Sphere of Influence (D < 5 Mpc) data taken from the cata¬ 
logue of lKarachentsev eta DEqqI. The Schechter function from Fig. ll5l is 
shown as the solid line and its extrapolation as the dotted line. The dashed 
line shows the contour of equal contribution to the overall luminosity den¬ 
sity. 


4.2 The surface brightness distribution 


As we f ound earlier th e Choloniewski function adopted by previous 
studies Jde Jong & Lacey .*2000!: fcross & DriveJll2002 ) provides a 
poor fit to the joint luminosity-surface brightness distribution. This 
is mainly due to the apparent broadening of the surface brightness 
distribution towards fainter luminosities. Fig. ll8l highlights this by 
showing the surface brightness distribution at progressively fainter 
absolute magnitude intervals along with the selection limits taken 
from Fig. 1121 We attempt to fit each of these distributions (Fig. 1181 
solid lines) with a Gaussian function defined as: 
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where, is the characteristic (peak) surface brightness and <r M e is 
the dispersion (analogous to those values in the Choloniewski func¬ 
tion). The errors reflect yTI-statistics from the original observed 
number of galaxies contributing to each histogram element. For el¬ 
ements with zero values (i.e., no galaxies detected) we adopt the er¬ 
ror appropriate for a detection of one galaxy so that the significance 
of these ‘null’ detections can be used to help constrain the fits. The 
number of degrees of freedom (listed in the final column of Tabled 
is the number of data values plus any zero values either side of the 
distribution within the selection boundaries minus the number of 
fitting parameters (three). Fig. 1191 a) shows the p e * (dark grey line) 
and 1 -cr M e values (shaded region) as a function of absolute magni¬ 
tude. These results are also tabulated in Tabled Fig. 1 191 a) shows 


a trend of an initially invariant luminosity-surface brightness rela¬ 
tion until Mb — 5 log h as —19 mag, and then a steady decline to 
the limits of detection at Mb — 5 log h m —14 mag. The disper¬ 
sion of the surface brightness distributions grows steadily broader 
as the luminosity decreases. The changes in both the luminosity- 
surface brightness relation and its dispersion explain the poor fit of 
the Choloniewski function which cannot accommodate such fea¬ 
tures. 

At this stage it is also worth noting that the Gaussian fits re¬ 
turn poor Xmin values at intermediate luminosities with even a hint 
of bimodality (at Mb — 5 log h ss —16 mag). This may imply that 
the low luminosity distribution consists of two overlapping pop¬ 
ulations: rotationally supported disks (lower surface brightnesses) 
and pressure supported spheroids (higher surface brightnesses). If 
borne out by more detailed structural studies then one might expect 
the rotationally supported disk systems to exhibit both a steeper 
M-p relation and a narrower surface brightness distribution than 
the overall population. 


4.2.1 Comparison to established structural trends 

Three luminosity-surface brightness trends (L-'S) are well known 
(see Table l2l : the Kormendy relatio n for spheroids iKormendxl 
1977 ): Freeman's Law for disks If reemat3 j 197 (|) : and the 
Virgo/Local Group dwarf rel ation | lmpey 1 Bothun & Mali nil 1 98^3: 
Mateo fl998l) . More recently Ide Jong & ulceTToPoll identified a 
relation for late-type disks similar to that seen for dwarf systems 
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Figure 18. Slices of the BBD along the surface brightness axis in intervals of constant luminosity as indicated (histogram with errorbars). Note that the data 
within shaded regions are not used in the determination of the Gaussian fits (shown as solid lines). 


in Virgo and the Local Group. The general trends for these estab¬ 
lished relations are overlain on Fig. ll9I a). One can speculate that 
the combination of these trends might well lead to the overall shape 
of the global L-E-distribution shown on Fig. ll9l a). as the spheroid- 
to-disk ratio typically declines with decreasing luminosity. Flence 
the global trend initially follows a Kormendy relation, switching to 
a more constant Freeman La w and finally into a decl ining dwarf 
relation. In comparison to the lde Tong & T.acevi 12000 ) sample for 
late-type disks we find a shallower luminosity-surface brightness 
relation and a broader surface brightness distribution. This may 
have two reasons: firstly, the dJL sample is /-band selected and 


our use of a global filter conversion may introduce this error; sec¬ 
ondly, the dJL sample is far more local and both inclination and 
dust-corrected. An obvious extension to the MGC analysis would 
be to also implement these corrections. 

Obviously, detailed structural analysis is required before a 
meaningful comparison to these ad hoc local studies can be made. 
However, we can highlight two key questions: (i) Are there two 
distinct disk populations, i.e. a giant disk population adhering to 
Freeman’s Law and a dwarf-disk population which does not, or 
just a single dwarf-disk relation? (ii) Is the dwarf surface brightness 
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distribution bimodal reflecting the two dominant classes of dwarf 
galaxies (dwarf ellipticals and dwarf irregulars)? 


4.2.2 Comparison to previous global L-E studies 

There are three significant studies of the global L -E distribution 
(see Table a 2dFGRS study by ICross et ail 1200 ll) . an SDSS 
study by IShe netal 120031) and a local galaxy SDSS study by 
iBlanton etalTl2005i ). 

In comparison to ICross et all 1200 ll) we find that the MGC 
data follow a similar L -E relation but show significantly broader 
surface brightness distributions. Fig. 1191 b) shows the comparison 
where we have taken the data from Table C2 of lCross et all1200 11 
and fitted Gaussian distributions in a similar manner as for the 
MGC data. The obvious discrepancy is in the values and this 
reflects a number of significant improvements. Firstly, the deeper 
MGC isophote coupled with the improved redshift completeness of 
MGCz (see Fig.0 enables the inclusion of lower surface brightness 
systems. As these preferentially occur at lower luminosity the effect 
acts to broaden the surface brightness distribution at fainter lumi¬ 
nosities. Secondly, the 2dFGRS sample was not seeing-corrected 
and hence the sizes of the high surface brightness galaxies were 
overestimated. Thirdly, the Cross et al. analysis derived the surface 
brightness measurements under the assumption that all galaxies 
were perfect exponential disks. These factors conspire to push the 
2dFGRS BBD into a narrower distribution. A fourth factor, hard to 
qualitatively assess, is the major revisions to the 2dFGRS photome- 
try which have taken place since Cross et alJ ( 200ill (see Cross et aH 


120041 for further details on the 2dFGRS recalibration history). 

In Fig.fTTlc) we overlay the recent SDSS results bv lShen et all 
( 2003 ). The Shen et al. study is ro.i-band selected with an imposed 
surface brightness cut at /r® [Circular] = 23 mag arcsec -2 . Shen et 
al. divide their sample into a number of sub-groupings and we elect 
to show the data from their Fig. 13 whic h separates galax ies accord¬ 
ing to concentration index, c (see Nakamura et all2003l for details). 
This divides the population into galaxies earlier or later then SO/a 
(with an 83 per cent success rate). To convert the Shen et al. r- 
band results to .Bmgc we derive the early and late-type (B M oc — r) 
colour in each of the r-band absolute magnitude bins specified by 
Shen et al. To do this we assume our eyeball E/S0 type corresponds 
to ‘early’ and our eyeball Sabcd/Irr to Tate’ (a detailed account 
of the eyeball morphologies and morphological luminosity func¬ 
tions is presented in Driver et al., in preparation). Hence each data 
point is adjusted independently. In addition we must also correct 
the SDSS effective surface brightness measurements from circu¬ 
lar to major axis measurements as discussed in Section 14.1.21 i.e. 
p e [MajorAxis] = p e [Circular] + 0.37. Fig. ll9r c) shows the com¬ 
parison between the Shen et al. data and the MGC and we see a 
uniform offset of 5p, e at 0.4 mag arcsec^ 2 . Despite extensive ex¬ 
perimentation with various data subsets and selection boundaries 
we have been unable to reproduce this result from the MGC data. 

Fig. 1 1 9l d) shows a comparison be tween the recent very low 
redshift SDSS r results ^Blanton et alJl2005al l and the MGC. We 
implement the same colour and surface brightness corrections as 
for the Shen et al. data. The ridge line of the very local SDSS 
data shows good agreement with the MGC data, albeit with a 
narrower distribution (i.e., lower o values). However this appears 
to be a further manifestation of the SDSS use of circular aper¬ 
tures. For example, if we derive the surface brightness distribution 
around the M*-point for the MGC data using circular and major 
axis derived surface brightness measures we find a peak offset of 
0.42 ± 0.03 (slightly larger than our estimated median offset of 


0.37) with UMajorAxis = 0.74 + 0.02 and a circular = 0.60 + 0.02. 
This latter val ue is now consistent with the recent SDSS estimate 
^Blanton et alll2005ah . Note that the dashed lines on Fig. ll9l indi- 
cated the SDSS extrapolated data used to correct the faint-end of 
the luminosity function. To the MGC limit these extrapolated dis¬ 
tributions appear to match our data well providing som e vindica¬ 
tion for the extrapolation process (see lBlanton et alJ2005aL Section 
l4T~3l and Fig .1771. 

_ The above two comparison b etween the MGC and the 

IShen et all 1200 3l) and iBlanton et all l2005al) studies must imply 
an inconsistency between the two SDSS results. This is perhaps 
more obvious if one bypasses the MGC data and directly compare 
Fig. 14 of Shen et al with Fig. 5 of Blanton et al. The peaks of 
the late-type/low sersic index SB distributions (both measured in 
r) are clearly distinct. Given the consistency with the 2dFGRS and 
the traditional structural studies discussed earlier we conclude that 
the problem most likely lies with the Shen et al. analysis. If we dis¬ 
regard the Shen et al study, we can conclude that the ridge line of 
the effective surface brightness distribution is well constrained and 
consistent between the MGC, 2dFGRS and the most recent SDSS 
results. However the MGC distribution is significantly broader than 
either the 2dFGRS or SDSS results which is due to improvements 
in the analysis (in comparison to the 2dFGRS) and the use of ellip¬ 
tical as opposed to circular aperture photometry (in comparison to 
the SDSS). 

4.3 Constraining the evolution of M* galaxies 

In Section 13.31 we described our adopted evolutionary parame¬ 
ter, j3, where the luminosity evolution follows the form: E(z ) = 
2.5 log[(l + z)~ a ], In Section l4.1.2l we saw that this evolution dif¬ 
fers from that adopted by the 2dFGRS and SDSS2 and discussed 
the possible impact. Fig. |20| (upper panels) highlights this depen¬ 
dency by showing the derived values for M* and a for various val¬ 
ues of f3. Taking this one step further we can attempt to constrain 
(3 directly by adopting the additional constraint that the recovered 
M* parameter should be invariant to redshift. We divide our sample 
into two redshift intervals: 0.013 < z < 0.1 and 0.1 < z < 0.18 
(with median redshifts of 0.08 and 0.13 respectively) and redeter¬ 
mine the collapsed LF distributions using our BBD analysis for 
regular intervals of (3. We hold a fixed to the value derived from 
the full redshift range for each (3 (i.e., the value indicated in the 
middle panel of Fig. HD and solve for M* and <j>*. It is neces¬ 
sary to hold a fixed because the high-z sample has too few low- 
luminosity galaxies to adequately constrain it. Fig. l20l (main panel) 
shows the resulting AM* versus j3 and we find the weak constraint: 
—2.0 < /3 < 1.25. While this is not particularly stringent, it is 
consistent with our adopted value of (3 = 0.75. This constraint on 
beta is margin ally inconsistent (2cr) with the evolution found by 
iBIanton et alJ (2 003a ) for the g-band ((3 at 2), but consistent with 
the evolution adopted by the 2dFGRS (/3 at 1). Perhaps more im¬ 
portantly this range of (3 implies a potential systematic uncertainty 
on the quoted M* (a) values of ±0.25 (±0.04). As noted at the 
start of Section^ it is also equally plausible that this error does not 
actually reflect evolution per se but could instead be interpreted as 
a redshift dependent photometric error or similar. 


5 CONCLUSIONS 

In this paper we have attempted to recover the space density and 
surface brightness distributions of galaxies, paying careful atten- 
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Figure 19. Four panels showing the 1-cr^e ranges (grey shading) and the fi e * values (dark grey line) from the Gaussian fits of Fig. flsl On all panels the 
MGC selection boundary is shown by the outer dotted line and represents a volume limit of 1800 h~ 3 Mpc 3 , (see Fig. 1 121 . (a) Overlain are the 1-cr r anges 
for the canonical spheroid (Kormendy relation; short dashed line), disk ( Freeman’s Law; so lid horizontal line) and late-type ide Jong & Lacevl|200(A long 
dashed line) res ults, (b ) A compa rison to the bj -band 2dFGRS results of|c ross et a D EooiS). data points and 1-cr ranges, (c) A comparison to the r-band 
SDSS results of IShen et alJ J2003h . Details of the colour correction are given in the text. Shen et al. divide their sample by concentration index into ‘early’ 
(c > 2.86, black symbols) and ‘late’ types (c < 2.86, grey symbols). At most magnitudes the late-types dominate the space density and so the most 
meaningful comp arison , band-sw itch not withstanding, is between the MGC data points and the SDSS late types, (d) A comparison between the global SDSS2 
distribution iBlanton et all2005alt transformed from r to -Bmgc as for (c). The dashed symbols indicate the region of extrapolation. 


tion to surface brightness selection biases. We have used the Mil¬ 
lennium Galaxy Catalogue which, although smaller in size, has sig¬ 
nificantly higher resolution, greater depth and higher redshift com¬ 
pleteness than either the 2dFGRS or SDSS datasets, hence probing 
to fainter intrinsic luminosities. In comparison to most earlier stud¬ 
ies this work includes the following enhancements. 


1. Imaging data which probes ~ 1-2 mag arcsec 2 deeper than the 
2dFGRS and SDSS surveys. 

2. Exceptionally high spectroscopic completeness (96 per cent). 

3. Individual K -corrections derived from uBgriz fits to spectral 
templates. 
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Figure 20. Main panel: constraining the evolutionary parameter /3 from the 
variation of M* between low and high redshift samples. Upper panels: the 
M* and a values for the entire dataset for different beta values. 

4. Seeing-corrected half-light radius measurements without as¬ 
sumption of profile shape. 

5. Modeling of the MGC surface brightness detection and reliabil¬ 
ity limits which includes the effects of the analysis software. 

6 . A weighting system which accounts for apparent magnitude and 
apparent surface brightness dependent redshift incompleteness. 

7. A joint luminosity-surface brightness step-wise maximum like¬ 
lihood method which incorporates the selection boundaries. 

5.1 The luminosity function 

Having incorporated these additional improvements to the clas¬ 
sical measurement of galaxy luminosity functions, we generally 
find close agreement with previous results. In particular, our col¬ 
lapsed luminosity distribution is marginally hig her in normalisa¬ 
tion a nd marginally flatte r than the 2dFGR5_[Norberg etak |2003l . 
ESP JZucca etal.l l 1997li and SPSS I iBIanton et all 1200 A after 
renormalisation bv iLiske' et alj l2003h results. It is significantly 
brighter and steeper than the more recent SDSS2 luminosity func¬ 
tion JBlanton et al.ll5003a ). We infer that this discrepancy arises 
because of a colour-selection bias in the SDSS. The MGC sur¬ 
vey extends to lower luminosity than these previous surveys, pro¬ 
viding more leverage and hence more reliability at the faint end. 
Our final Schechter function parameters are: <j>* = (0.0177 ± 
0.0015) h 3 Mpc -3 , Mg - 5 log h = (-19.60 ± 0.04) mag and 
a = —1.13±0.02 with ax 2 /v of 19.6/15 indicating arespectable 
fit. Overall we find that surface brightness selection effects do not 
play a significant role at bright luminosities. However, the combi¬ 
nation of the luminosity-surface brightness relation and the broad¬ 
ening of the surface brightness distribution is likely to lead to the 
loss of both low surface brightness and compact dwarf systems in 
contemporary surveys - hence the variation in a (cf. Fig. 1 1 5I & I 1 61 
and discussion in Section |4~H . At the very faint end we provide 
lower limits which suggest that the luminosity function must turn 
up at some point faintwards of Mg MGC — 5 log h « —16 mag. 


5.2 The low redshift bj luminosity density 

The final MGC Schechter function parameters, which include com¬ 
pensation for the effects of selection bias, imply a luminosity 
density of jbj = (1.99 ± 0.17) x 10 8 /iLq Mpc -3 (ad opting 
Mp, = +5.3 mag and Bmoc = bj — 0.131, cf. fL'iske et aH 
12003) . This is fully consistent with the MGC revise d values for 
the 2dFGRS, ESP and SDSS1 surveys (see Table 3 of ILiske et all 
2003). It is also worth noting that this value is also co nsistent with 
the o riginal 2dFGRS JNorberg et alJl2Qo3) and ESP JZuc caetalJ 
1 19971) values. The more recent SDSS2 result JBlanton et all2003al) 
finds jl=° = (1.54 ± 0.10) x 10 8 fiL 0 Mpc -3 and jg= 01 = 
(1.90 ± 0.10) x 10 8 /iXg Mpc -3 . Hence while the a = 0.1 
value is consistent the z = 0 value is not (3cr). The SDSS z = 0 
value is derived by adopting the evolutionary parameters derived by 
IBIanton et alJ J2003al) . If we use our evolution instead ( f3 = 0.75, 
see Section 13.3> the SDSS 2 = 0 luminosity density becomes 
jl =0 = (1.72 ± 0.10) x 10 8 fiZ /0 Mpc -3 which although closer 
is also inconsistent with our result. From Section 14. 1 .21 we con¬ 
cluded that a colour bias may have resulted in an underestimate of 
the SDSS g-band LF and hence to an underestimate of the SDSS 
luminosity densities. 

Based on the consistency between the MGC, 2dFGRS and ESP re¬ 
sults we can also find that while the more sophisticated BBD anal¬ 
ysis does recover some additional flux from the low surface bright¬ 
ness population, it does not constitute a significant change (consis¬ 
tent with the findings of lPrive J19991) . 

5.3 The joint luminosity-surface brightness relation 

Our recovered bivariate brightness distribution exhibits a clear 
‘ridge’ of data well separated from the selection boundaries down 
to Mg — 5 log h ss —16 mag at which point both very compact 
and highly extended dwarf systems will be missed. This indicates 
that the recovered Schechter function is generally robust to surface 
brightness selection effects for luminous galaxies but that the dwarf 
population(s) remain elusive. We attempted to fit a Choloniewski 
function JCholoniewskil [Ulis ). essentially a Schechter function 
with a Gaussian distribution in surface brightness, to the over¬ 
all distribution. It introduces three additional parameters: a char¬ 
acteristic surface brightness; a Gaussian dispersion; and a lumi¬ 
nosity surface brightness relation. We find that this function fails 
to provide a satisfactory fit due to the broadening of the surface 
brightness distribution with decreasing luminosity and a distinct 
change in slope of the luminosity-surface brightness relation at 
Mg — 5 log h fa —19 mag. 

5.4 The surface brightness distribution(s) 

We find that the surface brightness distribution of galaxies exhibits 
interesting behaviour, with a well bounded Gaussian-like distribu¬ 
tion at L* , qualitatively consistent with the combination of the Ko- 
rmendy relation for spheroids and a 2x broader Freeman’s Law 
for bright disks. At L* we find that the surface brightness dis¬ 
tribution is well described by a Gaussian with parameters: <j>* = 
(3.5+0.1) x 10 - 2 fi 3 Mpc -3 , /r e * = (21.90+0.01) mag arcsec -2 , 
cr M e = (0.77 ± 0.02) mag arcsec -2 with a \ 2 l v of 9.8/8. At 
fainter luminosities the normalisation increases, the characteristic 
surface brightness becomes fainter and the dispersion broadens un¬ 
til the distribution appears almost flat (or possibly bimodal). The 
Schechter function fit to the Gaussian normalisations (column 4 of 



































Table 3) provides a basic completeness correction and yields con¬ 
sistent Schechter function parameters to those quoted in Section 
im i.e. <j>* = (0.0184 ± 0.0008) h 3 Mpc~ 3 , Mg — 5 log ft = 
(-19.55 ± 0.08) mag and a = -1.11 ± 0.02 with a * 2 /v of 
26.0/12. 

In comparison to the 2dFGRS study of lCross et alJ 12001 1 we 
find a consistent luminosity-surface brightness relation but with a 
significantly broader dispersion in surface brightness. This is pre¬ 
dominantly due to the implementation of seeing corrected size, di¬ 
rect (i.e., profile-independent) half-light radius measurements and 
improved completene ss for extende d sy stems. A c omparison to the 
SDSS results of IShen et alJ 120031) and iBlanton et alJ l2005al) are 
complex because of the SDSS standard of circular aperture pho¬ 
tometry. However we find that our results can be reco nciled with 
those of Blanton *et alJ 12005 at) but not with those of IShen et alJ 
(2003). We also note that the SDSS use of circular apertures also 
leads to significantly narrower surface brightness distributions. 


5.5 Confrontation with theory 

Various studies iF all & Efstathiod 1 198(1 IDalcanton et afl ll997l : 

IMo et al 119981: Ide Jong & LacevfcOOOtlBell et all2003l) relate the 

Gaussian dispersion of the surface brightness distribution, or equiv¬ 
alently of the R e 4 distribution, to the dimensionless spin parameter 

of the dark matter halo, A = J\E\ 3G _1 M hak / ^Peebles 1969h . 
The basic concept is that the baryons in the disk reflect the an¬ 
gular momentum state (rotational or dispersive) of the dark mat¬ 
ter halo. While the luminosity relates to the mass, the disk scale- 
length relates to the specific angular momentum via the spin param¬ 
eter. Numerical simulations generally predi ct values for a\ n g, ~ 
0.56 ± 0.04 iBullock et aimOOlL see also ICole & Lacevlll99d) . 
This is most consistent with our observations at fainter luminosi¬ 
ties. However, at brighter luminosities we find a significantly nar¬ 
rower dispersion ai n g e = 0.36 ± 0.01. This may imply that 
brighter, more mature galaxies have evolved via mechanisms cur¬ 
rently beyond the standard numerical simulations (i.e., the baryons 
have decoupled from the dark matter halo and evolved away from 
the halo properties). Alternatively it may argue for significant non- 
directional baryon infall onto the halo, muddying any coupling be¬ 
tween the primordial baryons and the halo iMurali e" al]2002t) . Fur¬ 
ther possibilities might also include the presence of bulges in high 
luminosity systems whose properties depend more critically on a 
central supermassive black hole jSdk & ReeJ 19981) . 

One intriguing result of the simulations is the lack of vari¬ 
ation ofthe su rface brightness or size distribution with mass 
IVitvitska et alJ •2002). However, the more detailed studies which 
actually follow individual halos do imply a greate r variation in A 
for low-mass systems. For examDle. lVitvitska et all ( 2002 ) trace the 
history of three halos through the merging process. They clearly 
show that A initially varies significantly as the halos grow but that 
it becomes more resilient to mergers once the halos reach sufficient 
mass. Ultimately all three halos converge towards a common value 
at high mass. This seems plausible as one would expect smaller sys¬ 
tems to be more affected by mergers than more massive systems. 


4 R e is the half-light radius measured in kpc. It is related to surface bright¬ 
ness by /+ = M + 2.5 log(27rf? 2 )+ 36.57 and hence cr^e = 2.17cxi n g e . 
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5.6 Summary 

The surface brightness distribution of galaxies has been hotly de¬ 
bated for over thirty years since the advent of Freeman’s Law and 
the conjecture that it is wholly a selection bias iDisne 3G976). 
While in many minds th e debate is resolved through deep HI 
studies_(e ; g ; |Zwaaneyil2^JiH|)_and L^-a absorption line stud¬ 
ies IMo et all^99i1 ~ lcheri^^LanzettaJl2()()3l) . others remain loyal 
to Disney’s conjecture (e.g. O’Neil, Andreon & Culliandre 2003) 
and still argue for a flat surface brightness d istributio n at all lumi¬ 
nositi es l Spravb errv et al Jll997t lDalcantonllT99i3 : IDalcanton et alJ 
Il99'j) implying a significant amount of hidden mass. At least in 
O'Neil et al.’s case this stance derives from neglecting the clear 
luminosity-surface brightness relation seen in their data. The global 
surface brightness distribution may indeed be flat but the low sur¬ 
face brightness systems are predominantly dwarfs, providing little 
enhancement to the global luminosity, baryon or mass densities. We 
have attempted to address this debate through careful consideration 
of most selection issues. We probe to depths where the distribution 
is clearly bounded and where we robustly demonstrate we could 
reliably recover giant low surface brightness galaxies were they to 
exist in significant numbers. They do not. 

While we address the selection bias in the redshift survey as 
a function of apparent magnitude and apparent surface brightness 
other potential selection effects remain unexplored. For example, 
the redshift completeness may well be a function of spectral type, 
in which case the assumption that all galaxies of similar apparent 
magnitude and surface brightness are drawn from the same red¬ 
shift distribution will not be correct. In particular, inert low surface 
brightness galaxies which may have insufficient features to read¬ 
ily yield a redshift and may be missed entirely. We also do not 
at this stage decompose galaxies into bulge and disk components. 
This is clearly required by both the observations (which show two, 
or possibly three, distinct luminosity-surface brightness relations 
for spheroid and disk components) as well as by theory and sim¬ 
ulations which argue for distinct formation mechanisms for these 
components (merger and accretion respectively; e.g. Pierani, Mo- 
hayaee & de Freitas Pacheco 2004). We hence intend to extend our 
analysis in three ways: 100 per cent spectroscopic completeness; 
bulge-disk decompositions; and extension to longer wavelengths 
and in particular the near-IR which promises smoother, dust-free 
profiles. 
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Table 1. The breakdown of the survey contributions to the final MGC redshift sample. 


Survey 

No. of galaxies 
with Q z ^ 3 

No. of best 

redshifts 

Reference 

MGCz 

5054 

4944 

This work 

SDSS-DR1 

2909 

1528 

Abazaiian et al. ('2003') 

2dFGRS 

4129 

3152 

Colless et al. (2001) 

2QZ 

19 

13 

Croom et al. (2001) 

PF 

37 

28 

Francis et al. (2004) 

LSBG 

11 

2 

Imoev et al. (1996) 

NED 

1200 

29 

NED 


Table 2. Choloniewski (upper) and Schechter (lower) function parameters for the MGC galaxy distribution and comparable parameters from other surveys 
corrected to the Bmgc bandpass. MGC BBD refers to the empirical MGC BBD integrated over surface brightness. Values are for Qm — 0-3, =0.7 and 

H 0 = 100 km s -1 . 


Survey 

Filter 

r 

10~ 2 ft 3 Mpc —3 

M* 

^MGC 

mag 

a 

mag arcsec -2 

P 

mag arcsec -2 

Ref 

MGC (All types) 

^MGC 

2.13 + 0.32 

-19.37 + 0.15 

-0.99 + 0.14 

22.0 + 0.13 

0.33 + 0.15 

0.86 ± 0.09 

1 

Spheriods 

B 

— 

-19.70“ 

— 

20.03 

-1.57 

0.38 

2 

Disks 

B 

— 

-19.70“ 

— 

22.77 

0.0 

0.3 

3 

HDF (All types) 

Bab +0.13 

— 

-19.70“ 

— 

20.9 + 0.8 

0.67 + 0.09 

1 

4 

E/SOs (Total) 

R+ 1.4 

0.010 + 0.002 

-20.7 + 0.2 

1.06 + 0.26 

20.8 + 0.4 

0.92 + 0.09 

-0.62 + 0.08 

5 

E/SOs (Bulges) 

R+ 1.4 

0.008 ± 0.003 

-20.1 + 0.3 

1.30 + 0.50 

20.32 + 0.35 

0.84 + 0.09 

-0.55 + 0.10 

5 

Sd/Irrs (Total) 

Ic + 1-8 

0.332 + 0.07 

-19.7 + 0.1 

-0.93 + 0.10 

21.8 + 0.2 

0.51 ± 0.04 

0.61 + 0.04 

6 

Sd/Irrs(Disk) 

Ic + 1-8 

0.332 + 0.07 

-19.9 + 0.1 

-0.90 + 0.10 

21.7 + 0.2 

0.43 ± 0.05 

0.78 + 0.07 

6 

2dFGRS 

&2dF + 0.06 

2.06 + 0.09 

-19.67 + 0.02 

-1.05 + 0.04 

21.90 + 0.02 

0.28 + 0.02 

0.57 + 0.01 

7 

SDSS(Early) 

r 0 .i + 1.3 

— 

-19.70“ 

— 

19.4 

1.2 

0.6 - 1.2 

8 

SDSS(Late) 

r 0 .i + 0.9 

— 

-19.70“ 

— 

21.35 

0.54 

0.6 - 1.0 

8 

MGC BBD 

BmGC 

1.77 + 0.15 

-19.60 + 0.04 

-1.13 + 0.02 

— 

— 

— 

1 

ESP 

6j + 0.13 

1.65 + 0.3 

-19.59 + 0.08 

-1.22 + 0.07 

— 

— 

— 

9 

SDSS1 

g + 0.40 

2.06 + 0.23 

-19.64 + 0.04 

-1.26 + 0.05 

— 

— 

— 

10 

2dFGRS 

&2dF + 0.06 

1.61 + 0.06 

-19.60 + 0.008 

-1.21 + 0.01 

— 

— 

— 

11 

SDSS2o.o 

01 g - 0.09 

2.12 + 0.08 

-19.28 + 0.02 

-0.89 + 0.03 

— 

— 

— 

12 

SDSS2 0 .i 

01 g - 0.09 

2.18 + 0.08 

-19.48 + 0.02 

-0.89 + 0.03 

— 

— 

— 

12 


a Adopted value. ___ = ^^ _ _ ^_ __ 

References: (1) This work. (2)lKorm endy|{l977l ). (3) | Freemad |l970|). (4)|Drivei] j 19991) . (5)|de Jong et aU |2004i). (6)JdeJon^^Lace^j200C|). (7) ICross et alJ 
1200 lh . (8) IShen et alJ i20()3) . (9) IzuccaeT alJ 1 1997I) . (10) lBlantonetanT200lll . (11) lNorber^^in2002ll . (12) iBlantorTe^m^OCBal! " with and without 
evolution corrections. 
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Table 3. Gaussian fits to the surface brightness distribution in various intervals of absolute 
magnitude. Surface brightness is defined as the mean surface brightness within the half-light 
semi-major axis (i.e., using elliptical apertures). 


M Bmgc - 5 log h 

mag 

LL e * 

mag arcsec -2 

0> e 

mag arcsec -2 

r 

10 -3 /i 3 Mpc -3 

X 2 /" 

-21.25 

21.86 ± 0.10 

0.48 ±0.12 

0.08 ±0.02 

4.0/5 

-20.75 

21.72 ± 0.04 

0.58 ±0.03 

0.34 ±0.03 

14.0/10 

-20.25 

21.86 ± 0.02 

0.74 ± 0.01 

1.37 ±0.06 

5.67/10 

-19.75 

21.86 ± 0.01 

0.74 ± 0.01 

2.64 ± 0.09 

10.5/9 

-19.25 

21.94 ±0.01 

0.80 ±0.01 

4.33 ±0.10 

16.4/10 

-18.75 

22.10 ±0.02 

0.90 ± 0.01 

5.96 ±0.13 

30.5/11 

-18.25 

22.32 ± 0.03 

1.00 ±0.01 

7.29 ±0.22 

27.6/9 

-17.75 

22.58 ±0.03 

1.06 ±0.02 

8.86 ±0.27 

7.61/9 

-17.25 

22.86 ±0.07 

1.22 ±0.07 

9.69 ±0.65 

16.6/9 

-16.75 

23.00 ±0.09 

1.06 ±0.07 

9.48 ± 0.68 

22.7/8 

-16.25 

23.16 ±0.12 

1.20 ±0.10 

13.9 ±0.9 

16.0/7 

-15.75 

23.42 ±0.11 

0.86 ±0.12 

11.5 ±2.1 

4.12/7 

-15.25 

23.34 ±0.24 

1.32 ±0.20 

17.9 ± 3.0 

6.85/7 

-14.75 

24.26 ±0.43 

1.50 ±0.16 

19.8 ± 3.7 

4.31/6 

-14.25 

23.28 ±0.64 

1.14 ±0.40 

23.0 ±7.9 

1.38/3 


APPENDIX A: SIMULATIONS 

Galaxy sampling, detection and analysis selection effects are extremely complex (see e.g. Disney, Davies & Phillipps 1989) and depend on: 
galaxy profile; inclination; opacity; detector characteristics; survey depth and resolution; and the detection and analysis process, in particular 
background smoothing and removal. In this appendix we intend to crudely establish the empirical selection li mits of the MGC through 
simulations. The equations defining the redshift dependent selection limits are described in Section lT5l (see also lPrive Jl999l) and require a 
maximum and minimum size limit and a limiting surface brightness. In addition, we need a method to correct the observed half-light radii of 
compact galaxies for the effects of seeing. 

To establish the selection limits and the correction procedure we use IRAF’s ARTDATA package to simulate 45°-inclined optically thin 
exponential disks (this is consistent with our definition of surface brightness, see Section l±4l on a CCD frame with identical pixel size, sky 
noise and seeing as the MGC imaging data (i.e., 0.333 arcsec/pixel, ji LIM = 26 mag arcsec -2 and T = 1.27 arcsec). We produce simulation 
plates at 0.5 mag intervals from rn ir = 16 to 20 mag, each containing galaxies with half-light radii varying from r in = 0.1 to 32 arcsec (see 
Fis. lAlI . Compact galaxies were also simulated with de Vaucouleurs profiles but the results were the same as those for exponential disks. 

The MGC catalogue was derived using the SExtractor analysis software iBertin & Arnoutsil199(1) . The full details of this process are 
described in lLiske et al .. I S2003 1. We use SExtractor in an identical manner (e.g., sky background mesh 256 x 256 pixels with 7x7 filtering 
within each mesh) on the simulated data and also measure half-light radii in the same way as for the real data. 

Tables A1-A3 contain the results of our simulations as a function of m in and r in and these are also presented in Fig. lA2l 


Al Seeing correction 

From the top left panel of Fig. lA2l we can see the impact of the seeing, which was set to the median value of our survey. To crudely correct 
for this effect we make the resonable assumption that the observed half-light radius, r e , is given by the quadratic sum of the true half-light 
radius, r°, and some fraction, 77 , of the seeing FWHM, F: 

rl = (r°) 2 + t?r 2 . (Al) 

From the real MGC data we find that the median r e of stars is well described by 0.6 T (for a perfectly Gaussian PSF this would be 0.5 F) and 
hence one would expect 77 = 0.36. However, from the simulations we find rj = 0.32 and that equation <A 1 1 recovers the true half-light radii 
to an accuracy of better than 10 per cent for objects not classified as stars (see top right panel of Fig. lA2> . 


A2 Minimum size limit 

The most compact simulated galaxies are routinely classified as stars by the SExtractor star-galaxy separation and these objects are marked 
with an asterisk in Table Al. SExtractor's ability to distinguish a galaxy of given r e from stars cannot primarily depend on the absolute value 
of r e . It must rather depend on the galaxy’s distance, A, from the stellar locus in r e -space, where the scale of that distance is set by the width 
of the stellar r e distribution. From the MGC data we find the width of that distribution to be independent of seeing. Hence the minimum A 
at which a galaxy can still be reliably distinguished from stars cannot depend on seeing either. Recalling that the position of the stellar locus 
is given by 0.6 F we thus find: 


fmin — 0.6 r ± A MIN 


(A2) 
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Table Al. The recovered half-light radii (r 0 „i, uncorrected for seeing) from our simulations. Objects 
marked with a * were classified as stars by SExtractor. 


T*in 

arcsec 

16.0 

16.5 

17.0 

Simulated magnitude (mj n ) 
17.5 18.0 18.5 

19.0 

19.5 

20.0 

0.1 

0.63* 

0.63* 

0.63* 

0.64* 

0.64* 

0.64* 

0.67* 

0 .66* 

0 .68* 

0.2 

0.67* 

0.67* 

0.67* 

0.67* 

0.67* 

0.67* 

0.67* 

0 .68* 

0.70* 

0.4 

0.77* 

0.77* 

0.77* 

0.78* 

0.77* 

0.78* 

0.78* 

0.77* 

0.78* 

0.8 

1.07* 

1.07* 

1.07 

1.07 

1.07 

1.07 

1.07 

1.07 

1.07 

1.0 

1.27* 

1.27* 

1.25 

1.25 

1.27 

1.25 

1.22 

1.20 

1.20 

1.2 

1.40 

1.40 

1.40 

1.40 

1.40 

1.40 

1.38 

1.38 

1.38 

1.5 

1.67 

1.67 

1.67 

1.67 

1.67 

1.67 

1.67 

1.67 

1.65 

2 

2.13 

2.13 

2.13 

2.13 

2.08 

2.11 

2.08 

2.13 

2.16 

4 

4.00 

4.00 

4.00 

4.00 

4.00 

4.00 

3.89 

4.00 

4.00 

6 

5.55 

5.66 

5.66 

5.88 

5.77 

5.77 

5.77 

5.66 

5.44 

8 

7.99 

7.66 

8.10 

7.99 

7.88 

7.77 

7.44 

6.88 


10 

9.66 

9.66 

9.99 

10.21 

9.77 

9.54 

8.66 



12 

11.32 

11.32 

12.10 

12.21 

11.43 

10.99 




15 

13.65 

13.99 

14.77 

15.21 

14.32 

12.54 




20 

18.32 

17.98 

21.32 

18.87 

15.65 





22 

19.87 

19.42 

22.98 

20.42 






24 

20.42 

20.98 

25.31 

21.43 






26 

21.42 

22.53 

26.86 

22.76 






28 

23.09 

23.53 

28.97 

24.42 






30 

23.64 

24.64 

30.30 







32 

24.09 

25.20 

31.75 








To determine A M in we note from Table Al that in the simulations the confusion between stars and galaxies sets in at r M iN(T = 1.27) = 
1.07 arcsec, implying A„„ = 0.31 arcsec. The same value is found by fitting the actually observed minimum r e of galaxies in each MGC 
field as a function of T. 

Finally, combining equations <A 1 1 and <A2> we find the minimum seeing-corrected half-light radius: 
r°„ = \J 0.04 T 2 + 0.37T + 0.1 (A3) 


A3 Maximum size and low surface brightness limits 

Here we define these limits (on the seeing-corrected quantities) not as absolute detection limits but rather as the limits to which we can 
expect ‘reliable’ photometry. We arbitrarily define ‘reliable’ to mean magnitude and surface brightness estimates accurate to 0.1 mag and 
0.1 mag arcsec -2 respectively, thereby limiting the potential systematics due to isophotal errors and background subtraction to 0.1 mag and 
0.1 mag arcsec -2 . As most of our analysis uses 0.5 mag and 0.5 mag arcsec -2 -sized bins this ensures minimal impact of systematics on our 
final BBD. 

To determine r° AX and we proceed in the following way: using Tables A2 and A3 we find for a given m in the smallest r in for which ei¬ 
ther the recovered magnitude, m „ at , or the recovered (seeing-corrected) surface brightness differs by at least 0.1 mag or 0.1 mag arcsec -2 
from the respective input value. Given the resulting pair of ( m m , r,„) values we again use Tables A2 and A3 to find the corresponding (m onl , 
,) pair and plot this on Fig.|9](lower square points with errorbars). Finally, we find those values of and surface brightness which best fit 
these points. We find 

t"S*x = 15 arcsec (A4) 

and 

/x° IM = 25.25 mag arcsec -2 . (A5) 

Finally it is worth noting that some of the extreme low surface brightness galaxies were detected by SExtractor but they were broken up 
into numerous small objects. Rerunning SExtractor over the entire MGC with the deblender switched off may enable us to eventually probe 
fainter than the limits determined here, however a detailed study of the MGC using a cross-correlation method teoberts et all2004 failed to 
find any additional galaxies. 
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Figure Al. A simulated image showing 45°-inclined exponential profile galaxies all with magnitude 18 mag, but with half-light radii varying from 0.1 to 
32 arcsec (from left-to-right and bottom-to-top 0.1, 0.2, 0.4,0.8,1.0,1.2,1.5, 2, 4, 6, 8,10,12,15, 20, 22, 24, 26, 28, 30, 32 arcsec). X and Y units are in 
pixels (= 0.333 arcsec). Figure degraded see http://www.eso.org/~jliske/mgc/ for full pdf copy of this paper. 





26 


S.P. Driver el al. 


Table A2. The recovered magnitudes (m out ) from our simulations. 


Tj n 

arcsec 

16.0 

16.5 

17.0 

Simulated magnitude {m lTl ) 
17.5 18.0 18.5 

19.0 

19.5 

20.0 

0.1 

16.00 

16.50 

17.00 

17.50 

18.00 

18.50 

19.00 

19.50 

20.00 

0.2 

16.00 

16.50 

17.00 

17.50 

18.00 

18.33 

19.00 

19.50 

20.00 

0.4 

16.00 

16.50 

17.00 

17.50 

18.00 

18.49 

19.00 

19.52 

19.99 

0.8 

16.01 

16.51 

17.01 

17.51 

18.01 

18.50 

19.01 

19.51 

20.02 

1.0 

16.01 

16.51 

17.01 

17.51 

18.01 

18.51 

19.01 

19.53 

20.01 

1.2 

16.01 

16.51 

17.01 

17.51 

18.01 

18.52 

19.03 

19.53 

20.02 

1.5 

16.02 

16.52 

17.02 

17.52 

18.02 

18.52 

19.02 

19.51 

20.00 

2 

16.02 

16.52 

17.02 

17.52 

18.02 

18.52 

19.03 

19.54 

20.03 

4 

16.03 

16.53 

17.02 

17.53 

18.03 

18.55 

19.08 

19.58 

20.10 

6 

16.03 

16.54 

17.03 

17.52 

18.03 

18.55 

19.09 

19.64 

20.23 

8 

16.05 

16.56 

17.02 

17.50 

18.04 

18.58 

19.13 

19.73 


10 

16.06 

16.57 

17.00 

17.47 

18.02 

18.55 

19.16 



12 

16.09 

16.58 

17.00 

17.46 

18.03 

18.61 




15 

16.14 

16.59 

17.03 

17.44 

18.30 

18.68 




20 

16.11 

16.63 

16.89 

17.48 

18.17 





22 

16.13 

16.65 

16.88 

17.49 

12.13 





24 

16.23 

16.67 

16.89 

17.57 






26 

16.26 

16.70 

16.81 

17.51 






28 

16.27 

16.72 

16.91 

17.66 






30 

16.34 

16.79 

16.80 







32 

16.40 

16.78 

16.82 








Table A3. The input (^i„, left) and seeing-corrected output (^ ut , right) effective surface brightnesses from our simulations. Objects classified as stars are 
indicated by ***** and objects not detected are indicated by —. 


Tin 

arcsec 

16.0 

16.5 

17.0 

Simulated magnitude (mi n ) 

17.5 18.0 18.5 

19.0 

19.5 

20.0 

0.1 

13 00 ***** 

13 50 ***** 

14 00 ***** 

14 50 ***** 

15 00 ***** 

15 50 ***** 

16 00 ***** 

16 50 ***** 

17 00 ***** 

0.2 

14 50 ***** 

15 00 ***** 

15 50 ***** 

16 00 ***** 

16 50 ***** 

17 00 ***** 

17 50 ***** 

18 00 ***** 

18 50 ***** 

0.4 

16 01 ***** 

16 51 ***** 

17 01 ***** 

17 51 ***** 

18 01 ***** 

18 51 ***** 

19 01 ***** 

19 51 ***** 

20 01 ***** 

0.8 

17 51 ***** 

18 01 ***** 

18.51 18.50 

19.01 19.00 

19.51 19.50 

20.01 19.99 

20.51 20.50 

21.01 21.00 

21.51 21.51 

1.0 

18 00 ***** 

18 50 ***** 

19.00 19.06 

19.50 19.56 

20.00 20.11 

20.50 20.56 

21.00 20.98 

21.50 21.44 

22.00 21.92 

1.2 

18.39 18.40 

18.89 18.90 

19.39 19.40 

19.89 19.90 

20.39 20.40 

20.89 20.91 

21.39 21.38 

21.89 21.88 

22.39 22.37 

1.5 

18.88 18.91 

19.38 19.41 

19.88 19.91 

20.38 20.41 

20.88 20.91 

21.38 21.41 

21.88 21.91 

22.38 22.40 

22.88 22.85 

2 

19.50 19.53 

20.00 20.03 

20.50 20.53 

21.00 21.03 

21.50 21.47 

22.00 22.00 

22.50 22.48 

23.00 23.05 

23.50 23.57 

4 

21.01 21.00 

21.51 21.50 

22.01 21.99 

22.51 22.50 

23.01 23.00 

23.51 23.52 

24.01 23.99 

24.51 24.55 

25.01 25.07 

6 

21.89 21.73 

22.39 22.28 

22.89 22.77 

23.39 23.35 

23.89 23.81 

24.39 24.33 

24.89 24.87 

25.39 25.38 

25.89 25.88 

8 

22.51 22.55 

23.01 22.97 

23.51 23.55 

24.01 24.00 

24.51 24.51 

25.01 25.02 

25.51 25.47 

26.01 25.90 

26.51- 

10 

23.00 22.97 

23.50 23.48 

24.00 23.99 

24.50 24.51 

25.00 24.96 

25.50 25.44 

26.00 25.84 

26.50- 

27.00- 

12 

23.39 23.35 

23.89 23.84 

24.39 24.41 

24.89 24.89 

25.39 25.31 

25.89 25.81 

26.39- 

26.89- 

27.39- 

15 

23.88 23.81 

24.38 24.31 

24.88 24.87 

25.38 25.34 

25.88 26.07 

26.38 26.16 

26.88- 

27.38- 

27.88- 

20 

24.50 24.42 

25.00 24.90 

25.50 25.53 

26.00 25.85 

26.50 26.14 

27.00- 

27.50- 

28.00- 

28.50- 

22 

24.71 24.61 

25.21 25.09 

25.71 25.68 

26.21 26.03 

26.71- 

27.21- 

27.71- 

28.21- 

28.71- 

24 

24.90 24.77 

25.40 25.27 

25.90 25.90 

26.40 26.22 

26.90- 

27.40- 

27.90- 

28.40- 

28.90- 

26 

25.07 24.91 

25.57 25.46 

26.07 25.95 

26.57 26.29 

27.07- 

27.57- 

28.07- 

28.57- 

29.07- 

28 

25.23 25.08 

25.73 25.57 

26.23 26.21 

26.73 26.59 

27.23- 

27.73- 

28.23- 

28.73- 

29.23- 

30 

25.38 25.20 

25.88 25.74 

26.38 26.20 

26.88- 

27.38- 

27.88- 

28.38- 

28.88- 

29.38- 

32 

25.52 25.30 

26.02 25.78 

26.52 26.32 

27.02- 

27.52- 

28.02- 

28.52- 

29.02- 

29.52- 
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Figure A2. A graphical representation of Tables Al, A2 and A3. The top panels show the ratio of the output and input half-light radii both with (right) 
and without (left) seeing correction. These panels demonstrate the effectiveness of the seeing correction. The dashed lines show the minimum size limit 
v MIN = 0.6T + 0.31 arcsec. The middle panels show the difference between the output and input magnitudes as a function of the seeing-corrected output 
half-light radius (rjj,,, right) and as a function of the seeing-corrected output effective surface brightness (/ij ut , left). Finally, the bottom panels show difference 
between the seeing-corrected output and input surface brightnesses as a function of the same two parameters. The dashed lines in the left panels show the low 
surface brightness limit /zp IM = 25.25 mag arcsec -2 , while in the right panels they mark the maximum size limit r® AX = 15 arcsec. In all panels we show 
objects identified as stars by SExtractor as open star symbols, objects beyond rJJ AX as crosses and objects beyond /zp IM as open squares. Solid symbols represent 
galaxies within the selection limits and for these we recover magnitudes and surface brightnesses that lie consistently within ±0.1 mag and ±0.1 mag arcsec -2 
of their input values. 
















